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1.0 INTRODUCTION 


1 . 1 General 

Requirements for better performance and longer life have pushed engine 
designs to lighter weight systems, higher reliability, and increased 
pressures and environments. Temperatures, external and internal fluid flow 
noise, and mechanical vibration levels have increased markedly and have been 
shown to limit the hardware designs. Advanced engine concepts and designs 
are different enough that the loads cannot be simply scaled from other 
engines. 

The use of engine cycles such as staged combustion on the SSME result in 
engine operating pressures in the 3000 to 7000 psi regime. High performance 
turbomachinery operate in the 30,000 to 100,000 RPM regime. These 
operational requirements result in complex high energy loading throughout 
the engine. The difficulty in installation, cost, and the potential for 
destroying an engine has severely limited the required instrumentation and 
measurements to adequately define loads of key components such as turbine 
blades. Also, accurate analytical methodologies for defining internal 
flow-related loads are just emerging for problems typically found in rocket 
engines. The difficulty of obtaining measured data and verified analysis 
methodologies has led to the probabilistic load definition approach of this 
contract. 

Current loads analyses methodologies are driven by their usage in 

deterministic analysis methods. This includes strength and fatigue analysis 
as well as mechanical vibration. The deterministic solution typically uses 
an upper bound approach where maximum loads and minimum properties are 
used. For critical hardware, a separate sensitivity study is often made to 
determine more nominal operation and which loads and their variation govern 
the hardware design. 
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The Composite Loads Spectra Contract (CLS) and the associated Probabilistic 
Structural Analysis Method (PSAM) contract from Lewis Research Center are 
developing an integrated probabilistic approach to the structural problem. 
The probabilistic loads approach has the ability to more technically 
quantify knowledge relative to the loads. The use of mean values and 
distribution about this central value rather the maximum or enveloped loads 
can add greatly to the understanding of normal engine operation and still 
furnish as good or better knowledge of maximum conditions. 


The present techniques often results in manufacturing of components that in 
many cases greatly exceed design requirements, but there is no way of 
assessing this margin for extending the useful life margin. Thus, to 
formulate more effective designs, it is necessary that the loads on the 


components of rocket engines be derived so that they can be applied by 
probabilistic analysis methods such as PSAM to end up with results that are 
quantifiable to more accurately reflect the true risk. The SSME engine is 
currently undergoing a failure modes affect analysis. The assessment would 
be much easier to perform if a probabilistic analysis and associated risk 


assessment were available. 


This project will provide methods to combine technologies of analytical 
(deterministic) loads and probabilistic modeling. Since these methods will 
be developed from a generic approach, they will be applicable to current or 
advanced liquid rocket engine designs. 



The objective of this program is to develop generic load models with 
multiple levels of progressive sophistication to simulate the composite 
(combined) load spectra that are induced in space propulsion system 
components, representative of Space Shuttle Main Engines (SSME), such as 
transfer ducts, turbine blades, and liquid oxygen (LOX) posts and system 
ducting. The approach will consist of using state-of-the-art probabilistic 
methods to describe the individual loading conditions and combinations of 
these loading conditions to synthesize the composite load spectra simulation. 
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The methodology required to combine the various individual load simulation 
models (hot-gas, dynamic, vibrations, instantaneous condition, centrifugal 
field, etc.) into composite load spectra simulation models will be developed 
under this program. Results obtained from these models will be compared 
with available numerical results, with the loads induced by the individual 
load simulation models, and with available structural analysis results from 
individual analyses and tests. These theories developed will be further 
validated with respect to level of sophistication and relative to predictive 
reliability and attendant level of confidence. 

A computer code incorporating the various individual and composite load 
spectra models is being developed to construct the specific load model 
desired. The approach is to develop incremental versions of the code. Each 
code version will add sophistication to the component probabilistic load 
definition and the decision making processes, as well as installing a new 
set of loads for an additional component. This allows for ongoing 
evaluation and usage of the system by both Rocketdyne and NASA. 
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2.0 SUMMARY 


2.1 General 

The development of probabilistic generic load models is a 3-year base 
program and a 2-year option program. Rocketdyne is responsible for the 
overall project. Battel le Columbus Laboratories is the major subcontractor 
for developing probabilistic load models and furnishing technical expertise 
in probabilistic modeling in general. 

The effort is divided into three tasks: the probabilistic model theory, 

code development, and code validation and verification. An initial survey 
effort was made to review available L0X/LH 2 data on the components under 
study and appropriate probabilistic load methodologies for use in this 
contract. Four rocket engine components, LOX posts, transfer ducts, turbine 
blades, and an engine system duct are being used as example components for 
the loads development. Examples of these components are shown in Figures 
1-3. Figure 1 is a cross section of the SSME powerhead showing the LOX 
posts in the 3 combustors and the transfer ducts between the powerhead 

components as well as the standard instrumentation that is used for 
monitoring the engine. Figure 2 shows a typical turbopump with its two sets 
of turbine blades. Figure 3 is an overall SSME powerhead view where the 
system ducts are depicted. Of specific interest is the high pressure 

oxidizer turbopump discharge (HPOTPD) duct. 

Simply stated, the goal of the composite load spectra project is to provide 
a tool to generate probabilistic based composite loads of a rocket engine 
design. These loads can be used to improve aspects of current deterministic 
analysis approaches or as input to a probabilistic analysis method such as 
PSAM. In the first year, an initial code was developed that had the 

essential features of the planned expert system and probabilistic loads. 
This code was limited in scope to steady state turbine blade load components. 
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Individual 
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Components 


turbine blades were identified at the start of the project as specific 

components for load development. The fourth component chosen for this 

project was an engine system duct, the HPOTPD duct. The oxidizer ducting 

system on the SSME has experienced a series of problems related to flow 
vibration that were unexpected. High energy flow vibration environments and 
their application to hardware analysis have not been well developed to 
date. By choosing this component, additional load definitions will be 
developed to aid in understanding and minimizing potential problems in 

future rocket engine designs. 
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2-2 Probabilistic Load?; Development 

One of the goals of the program is to be able to address generic engines 
that may include different mission profiles or incorporate design changes. 
This requires that a robust and general probabilistic approach be adopted 

for inclusion in the expert system model. During the first year of the 
program, a survey was conducted to select these probabilistic models and the 
initial programming, debugging and shake-down analyses were performed. The 
second year of the program has been oriented towards refining the 
methodology, developing a database that can be used by both the 
probabilistic methodology as well as the expert system, including different 
functional forms for the load description, model verification and 
validation, and the generalization of the computer program system. 

The probabilistic model has included three probabilistic methods: 1) a 

second statistical moment propagation method which assumes that all of the 
load variables and engine parameters are normally distributed, 2) a 

discrete probability method (RASCAL), and 3) Monte Carlo analysis. The 

moment propagation method, referred to as the Quick Look Model (QLM) 

provides a fast, efficient method for determining the composite load 
distribution if the basic distribution of variables are not severely 

skewed. The RASCAL method is a discrete method capable of handling standard 
distributional forms, e.g. normal * lognormal, Weibull, and so on, 

non-standard forms such as bi-modal, and provides a range of levels for 
accuracy. This method can also be used to perform importance sampling which 
can be used to examine regions of concern for the composite load even though 
such values are rare. Finally, Monte Carlo analysis is available so that 

classical confidence limits can be obtained to assess the accuracy of the 
composite load prediction. 

All phases of the mission history profile are addressable by the 

probabilistic load model.' Currently, each portion of the mission history is 
defined as transient, quasi-steady, or steady state phases. The transient 
phase is characterized by rapid changes in the amplitude of the individual 

loads and engine parameters. 
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The rapid changes allow the program to ignore small oscillations about the 
much larger load fluctuations. The uncertainty in the load enters from the 
variability in the peak load value and its time of occurrence. The 
quasi-steady phase is that portion of the mission where the nominal value of 
the load is slowly changing, and thus can be approximated by "staircase- 
type steady state steps. The steady state region is where the nominal 

values of all of the individual and composite loads are constants. Unlike 
the transient phase, both the quasi-steady and steady state phase do have 
fluctuations superimposed upon the nominal behavior. Additionally, each of 
these phases can have "spike" values superimposed which represent the 

occurrence of rare events. 


The linking of these different mission phases has been completed. It has 

been demonstrated that for the cases where data have been available that a 
continuous, nominal behavior is achieved. In addition, the Predic .d 
variability and the measured variability are well within acceptable limits. 
Therefore, the extension of the model has proceeded to engines and mission 
definitions for which little or no data exist. 


For the SSME engine, expert opinion data was obtained by Rocketdyne for 
those loads and engine parameters used in the model for which measurements 
were unavailable. Test runs of the probabilistic model with these data 
included were made and compared to measured composite load data. The 

results indicated that the variability in composite load type data was 
adequately predicted by the model. Some differences in the predictions and 
measurements have been found and will be further looked at as part of the 
validation phase of the work. Late in the year, some analyses were begun 
which examines other engine types. The approach for determining the mean 
and variability in the individual load parameters for engines for which no 
data exist is to scale them using the engine design parameters and table 
look-up values developed during the second year of effort. Further 
development and validation of these analyses will be performed during the 

third year of the program. 
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Documentation of the code, ANLOAD, has continued throughout the program. 
Periodically, new versions of the program are sent to Rocketdyne for 
incorporation into the expert code system. The code work to date has 
primarily addressed loads that are dependent on the overall engine 

performance and are directly relatable to the engine model and duty cycle. 
The next phase of the loads development will address the remaining loads, 

e.g. fluid and mechanical vibration environment, sideloads and shock, pops 
and chugs and debris loads. Specific modeling for each of the components 
will also be completed. Additional specific modeling of the four components 
is also required. Representative loads for each of the four components in 
the study will be used for the validation and verification of the code. 

2 . 3 Load Expert System 

The probabilistic loads model is implemented as part of an expert system. 
The expert system is a tool to generate and analyze composite loads of a 

rocket engine design and to supply these loads for use in either 

deterministic or probabilistic FE computer codes for performing structural 
analysis of engine components. The statistical information used in the 
expert system primary basis is SSME test results, but expert opinion and 
other available engine data are used when appropriate. The approach is to 
develop the knowledge base of an individual load formulation on a reasonable 
physical basis in as generic a sense as possible. Engine statistical data 
are part of the knowledge base and used where appropriate. 

A knowledge-based system has the facility of building up a large domain 
knowledge base and maintaining a large amount of data. It has the 
capability to perform logical deduction and inferences and thus it can help 
users to make decisions and to solve problems. These characteristics allow 
one to build an expert system to simulate and perform the process of 
problem-solving by an expert in a particular problem domain. 

The functions of this knowledge-based system are to manage the database, 
provide expert knowledge in generic probability loadings for rocket engine. 
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A FORTRAN based non-proprietary knowledge system development tool that can 
satisfy all the needs of this project is not available. Therefore, it was 
decided early that the knowledge system will be built to suit the need of 
this project. 

A simple philosophy discovered by pioneering workers in the field is that 
the power of a knowledge base system is in its capability to have a vast 
amount of domain knowledge and not necessary to have a complex inferencing 
engine. Following this philosophy, the load expert system LDEXPT was built 
with a simple inference system. The expert system uses the ANLOAD module to 
perform probabilistic modeling and statistical analysis. To make knowledge 
representation more efficient for the load expert system, a database system 
was implemented. This database system facilitates the communication between 
the expert system and the knowledge base, helps to maintain data integrity 
and avoid data redundancy. The load expert system LDEXPT version 2.0 has 
all three elements in place. Its knowledge base has load information for 
SSME type engines, knowledge about the influence coefficient method based on 
engine performance analysis and initially the turbine blade load information 
and scaling model calculation. 

The load expert system is a rule-based expert system. The inferences are 
carried out with rules. In the load expert system, the rules are 
modularized. Each module is designed to solve a particular problem or to 
perform a task. The load expert system LDEXPT version 2.0 has rule modules 
to calculate turbine blade loads using scaling model and generate engine 
dependent loads (e.g. HPFTP discharge pressure) using influence coefficient 
method. The rules designed so far are mostly related to process control and 
information retrieval. In the next development, rules to generate 
probability models for a complicated composite load spectra will be designed 
which will require more use of artificial intelligence. 

The load expert system now has knowledge of the turbine blade loads for 
generating steady state and quasi-steady state load spectra. Additional 
load data on pressures and temperatures are ready for adding to the rules. 
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The transient loads, pops and chugs and vibration loads, etc. are being 
developed and will be implemented as soon as the model development is 
complete. Knowledge on the transfer duct has been collected and rules for 
transfer duct load calculations can now be developed. The other two 
components loads will follow. 

The basic expert system components of the load expert system LDEXPT: the 
expert system driver, the database system, the FORTRAN data management 
system and the basic probabilistic modeling and statistics tool box are all 
in place, that is the main tasks of system development phase are complete in 
version 2.0 of the code. The next main task is the further development of 
applications of the expert system to the composite load spectra project. 


13 



3.0 ENGINE LOADS 


3.1 Background 

The individual loads applicable to the four components in this project are 
summarized in Table 1. These loads cover a major portion of the loading 
throughout a rocket engine and are an excellent representative set to 
develop into an engine loads expert system. Where applicable, the 

individual loads are modeled for the entire duty cycle. 

The loads are essentially self-generated or induced loads except for steady 
state g-forces and gimoaiing requirements during flight. This allows the 

engines to be readily separated from the vehicle loads analysis as a 
subsystem with specific requirements. 

The vehicle design can be divided into conceptual, preliminary, detail and 
design verification phases. This is followed by flight support and possibly 
uprating and problem resolutions. During the conceptual and preliminary 
design phases of a vehicle, major decisions are reached that spawn 
requirements for engine design. Vehicle requirements often are related to 
load alleviation or preventative measures and performance requirements to 
optimize vehicle design with engine design. Examples are. 1) controlled 
thrust rise rate, 2) in flight load alleviation, 3) cutoff impulse 
requirements, 4) engine inlet operating pressures and temperatures, and 
engine gimbal angle and rate requirements. A description of the approach to 
deriving loads design criteria for the space shuttle and it s payload is 
given in Reference 1. The vehicle system requirements reduce to a set of 
engine loads and system requirements. Reference 2, that define limits and 
engine duty cycles that end up defining a part of the engine individual and 
composite loads. (Note: most of the examples in the discussion herein 

presented are related to the SSME, but it is appropriate to rocket engines 
from a generic standpoint.) 
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The basic engine duty cycle is controlled by engine thrust buildup limits, 
Figure 4, engine thrust decay limits, Figure 5, and overall flight 
requirements such as maximum operational power, throttling during maximum 
dynamic loads and throttling near the end of flight to maintain a maximum 
g-load limit, Figure 6. These and other engine requirements are used to 
develop engine configurations and models. The engine models furnish 
inter-related deterministic loads — pressures, temperatures, vibration 
levels, etc., for major components such as inlet and outlet conditions for 
transfer ducts, preburners, injectors, and turbopumps, see Figure 7. These 
interface loads are used with deterministic models to evaluate loads on 
individual components like turbine blades transfer ducts, LOX posts, etc. 
For instance, the steady-state loads are used by the hydrodynamics 
specialists to determine loads across each turbine stage or blade. The heat 
transfer specialists use information from the same model results to 
determine blade temperatures. The dynamics experts use the model results to 
determine turbine blade dynamics. The structural and analysis experts use 
the model information and the input loads from the other experts to develop 
the total load and structural analysis. 

Deterministic models of varying complexity are used in all analysis 
efforts. The steady-state engine simulation model can furnish discrete 
values at an operating point. The influence coefficients relate one or 
several engine parameters versus . other parameters. Somewhat similar 
information can be determined from the engine transient simulation model. 

Simulation models are formulated using generic engine process descriptions 
and constitutive equations and detailed tabulation of the propellant 
physical properties. The description of the basic processes of the system 
simulation involves all static and dynamic formulations (where applicable) 
that are considered of importance in accurately representing the overall 
behavior of the engine during start, mainstage control, and cutoff. The 
validity and veracity of these process descriptions in terms of their 
ability to describe the overall system behavior have been proven by 
correlation of simulation results with engine test results from previously 
developed rocket engines. 
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Figure 6. SSME Typical Flight Profile 


With appropriately defined changes in the coefficients of the process 
descriptions, any new or modified engine component can be modeled into the 
simulation. Thus, the analytical description of even an entirely new engine 
system can be formulated and used with a confidence level that is based on 
previous proven performance of the analytic basis. 

The engine performance model is a complex code not readily usable for the 
CLS effort. But engine influence coefficients are typically developed for 
rocket engines based on the performance model and are a practical method to 
develop a subset of the loads. Using an influence coefficient approach for 
the general operating conditions allows generic loads development across 
significantly different engine cycles. The three production L0X/LH 2 
flight engines developed by NASA have had different engine cycles: the RL10 
has an expander cycle, 0-2 had a gas generator cycle and SSME has a staged 
combustion cycle. 
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Figure 7. Interrelation of SSME Analysis Models 


In support of this project, SSME influence coefficients have been extended 
to relate key engine variables to additional turbine and hot gas system 
related parameters. These coefficients are applicable to all four SSME 
components addressed by this project. 

For the CLS work, the generic engine cycle is divided into start, cutoff, 
quasi-steady state and steady state operation. This operational mode will 
be discussed first. 

3.2 Steady State and Quasi Steady State Operation 

Except for transient conditions, nominal generic duty cycle loading can be 
described by a relatively few independent parameters — the power level 
variation and other engine direct variables, such as inlet pressures and 
temperatures. Using these independent parameters with the applicable 
influence coefficient, nominal operation conditions are readily determined 
at component interfaces throughout the duty cycle. 
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Variations about the nominal condition for a specific load or parameter is 
approached in two separate methods — one based on estimated engine random 
variation, and the other based on measured engine data. The estimated 
engine variation for selected independent variables (41 variables in the 
case of the SSME) can be used with the engine performance model to determine 
variations throughout the engine. These variations were developed from 
consultations with the individual experts on specific hardware and covers 
geometric, performance, etc. conditions that can effect the engine operation. 

Typically, an engine performance data slice is obtained for each engine test 
and flight after the engine operation is stabilized and at a consistent time 
period, e.g. 190 to 200 seconds after start for the SSME. This information 
is utiiizea to calculate basic engine performance and the engine, to engine 
and test to test variation of engine operation. Similar data is available 
for other engines such as the 3-2, Atlas, F-l , etc. 

The purpose of developing these variations is to furnish operating ranges 
for engine performance parameters and for use in validating the engine 
variations used in the model. These random variations are added to the 
predicted performance effects of direct independent variation allowed by 
specifications to determine parameter maximum and minimum expected values. 
The same percentage random variables are used throughout the thrust limits 
of the engine. The inherent assumption in this calculation of perturbed 
engine operations is that the 41 variables are independent random variables 
with normal distributions. 

The basic perturbation technique of the engine model used for this analysis 
is similar to that used for calculating the influence coefficients that are 
being used for the Composite Loads Spectra (CLS) work. For the influence 
coefficients, a matrix of the individual effects are maintained to allow for 
direct determination of individual variable changes. 

The measured engine data is only a small subset of calculated engine 
variables, but can be used to substantiate the calculated variations. 


19 



The 41 variables have counterparts, in general, to the 26 variables used in 
the engine influence coefficients. They are not identical since they are 
defined for different purposes. The 41 random variables as mentioned above 
are to cover all engine to engine and test to test variations for use in 
component design. The influence coefficients were developed for the 

customer's use in accounting for flight performance variations of a specific 
engine. 

This information is the best data currently available for use on the CLS 
contract. Currently, there is ongoing work to develop a set of 2 sigma 
variations of measured parameters (two standard duration bounds) 
specifically based on the engine test database, but this will not be 

available for several months. 

The direct independent variations include: propellant inlet temperatures 

and pressures, line resistance changes due to gimbaling, and tank 
repressurization flow settings. These maxima and minima define the 
operational limits used for engine component design. The engine ICD 
(Interface Control Document), e.g. Ref. 2, defines the required operational 
bounds of inlet pressures and temperatures that the engine must operate 
within. The gimbaling limits are also furnished in the ICD that were used 
to develop in-line resistance calculation input set. The effect of the 
direct variables are obtained by developing maximum non-compatible load 
variations based on the operational bounds -- the corners of the operational 
boxes. As mentioned above, these maximum non-compatible load variations are 
used as additions to the random variable perturbations to determine a 
maximum and minimum engine balance condition. Surge and transient effects 
are added as additional perturbation effects. 

The random variable perturbation information is consistent with the CLS 
approach and has been used in the probabilistic load model development for 
AN LOAD. The direct independent variations are duty cycle load parameters 
along with power level that defines part of the component loads. The direct 
variations have nominal values and perturbations based on engine data — see 
discussion in section "Comparison of ANLOAD Predictions with Expert Opinion". 
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The transient phase of the load definition is based on a combination of 
vehicle requirements, engine simulation models and engine test results. 
Typical vehicle requirements were discussed earlier — start and cutoff 
transient envelope 0 specified to minimize vehicle loads. Additional 
requirements like rates of power level changes during throttling and 
associated up-thrusts are additional requirements that size control system 
variables. These system requirements indirectly control some of the nominal 
loads on components during transient operation. Thrust control drives pump 
speeds, torques, and system pressures and temperatures. 

Various transient models are employed according to the type and range of the 
system dynamics under study. The analog model is generally used for 
surveying system characteristics, tradeoff and optimization of the control 
system, and in evaluating the large number of system changes typical of the 
early phase of engine design. Hybrid simulation is used to study digital 
control operation with the analog model. The hybrid computer thus simulates 
the role of an engine interfaced with a digital control system. The digital 
model, which most accurately represents system behavior, is used for 
simulation studies where maximum accuracy of results is needed, or where 
wide-range nonlinear operating conditions exceed the normal capabilities of 
analog simulation. 

The dynamic simulation models and steady state performance models describe 
the same processes, but the performance models stress accuracy of steady 
state operation parameters, whereas the dynamic simulation models have to 
consider the overall system behavior throughout the duty cycle. From a 
loads definition standpoint on this project, the SSME dynamic digital 
simulation model results are used for transient conditions below 657. power 
level and performance model results above this thrust level. 

The performance and dynamic simulation models are deterministic solutions. 
The transient solutions are essentially a nominal operation description of 
the engine operation. Engine to engine and test to test variation, as well 
as certain high frequency transient conditions or non-uniform flows and 
temperatures, are not adequately modeled for local load definition. Actual 
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engine measurements typically are not at the point of interest to define the 
transient and steady state operation conditions. Special instrumented 
components aid in this load definition, e.g. instrumented turbines. They 
may be close to the component in question, have better response 
characteristics, but usually survive or are utilized for a very limited 
number of tests. The definition of individual component load distributions 
require a combination of: 1) expert knowledge from previous engines and 

testing, special measurements, standard measurements and simulation models 
specifically formulated to calculate an engine test operation using measured 
conditions. 

The hot gas transient load distributions have been based on the SSME HPFTP 
hot gas side of the engine. A simulation model was constructed of the fuel 
side where a combination of engine measurements and the simulation model 
were used to define hot gas system and fuel turbine start and cutoff 
transients. Measured parameters included pump speed, turbine discharge 
temperature, and pump delta pressures. Turbine torque was developed from 
pump head and torque curves with corrections for initial torque of the 
turbopump. The transient temperature ignition spikes were based on a 
correlation of instrumented turbine temperature measurements and the 
standard temperature discharge bulb measurement from measured data. After 
the temperature spikes subside, the turbine inlet temperature was based on 
measured discharge temperature corrected for the heat loss across the 
turbine due to the work energy extracted. Using this methodology, a series 
of engine tests were processed. The tests selected covered the expected 
bounds of the high pressure fuel turbine system operation. The results of 
this study was used in developing the transient model using this turbine as 
an example. A test by test tabulation of HPFTP turbine temperatures of all 
SSME hot fire tests and flights was also used in developing a statistical 
database of expected turbine temperature variation. The database included 
start transient temperature spikes as measured by the turbine discharge 
temperature. The magnitude of the inlet temperature spikes for the database 
tests was calculated using the same correlation procedure described above. 
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3.3 Structural Dynamic Excitation 


The structural dynamic excitations used for rocket engine components are 
typically measured responses to combustion processes, turbomachinery 

generated loads or aerodynamic internal flows in ducts or nozzles. 

Accelerometers are located on external structure of major components that 
generate these loads such as the location(s) shown on the SSME main injector 
(the LOX dome interpropellant plate or the gimbal bearing flange connection 
and turbopumps, Figure 1). The accelerometers are standard measurements on 
test firings and engine flights. General vibration environments and engine 
redline limits are defined using the standard flight instrumentation. 
Additional accelerometer measurements are made for developing specific 
vibiation environments to be used on individual components. The measured 
responses are used as dynamic base input accelerations for individual 
components like a LOX post or transfer duct, or as input accelerations to an 
injector assembly model with the entire set of LOX posts, interpropellant 
plate and LOX dome, etc. The current state of the art is to use the 
response as an input rather than transform the responses back to the actual 
load functions. Accelerometer data is measured in one, two, or three 
mutually orthogonal directions and furnish local magnitude and frequency. 
This data is insufficient to identify the various mode shape, so simplifying 
assumptions are typically made. These include independent assessment of 
vibrations by load direction and the assumption that there is no correlation 
between accelerometers. 

The generic vibration mission-history-profile is complex since it can be 
made up of several different load components whose significance is variable 
and dependent on engine and component parameters (J-2, SSME, OTV, etc.). 
Pictorially, a typical vibration response mission-history-profile is shown 
in Figure 8. The loads can be categorized as: 

1 . Transient Loads 

a. Random pops (high frequency shock) - local combustion detonations 
during start and cutoff and up to minutes after cutoff. Pops can 
occur infrequently during the initial steady state condition. 
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Figure 8. Pictoral Representation of Generic Vibration Response 


b. Engine side load reactions (low frequency oscillations) - overall 
structural loading from the nozzle exhaust plume separation that is 
reacted by the primary load path through the engine structure and 
gimbal bearing and gimbal actuators. 

c. Nominal vibration - energy that builds up with the magnitude of the 
combustion-related engine power level and flows in turbopumps. The 
vibration level varies as the engine power level is changed 
throughout the duty cycle. 

2. Steady State Operation Loads 

a. Nominal random vibration - combustion and turbomachinery related 

mainly from the load generator nearest the accelerometer, but 
potentially from other load generators on the engine. 
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b. Steady state sinusoidal vibration - significant discrete sinusoidal 
vibrations are measured at multiples of pump speeds on turbopump, 
preburner and main injector accelerometers. 

The extensive engine test measurements have been taken with the standard 
accelerometers on virtually every engine test. The signals are processed with 
AMS/RMS, ISOPLOTS and STATOS records. Vibration levels and pops are tracked 
on a test-by-test basis. 

Zonal shock and vibration criteria are defined for the entire engine. The 
methodology currently used for defining the loads envelopes the maximum 
responses from at least three tests each on two engines at the power level 
within a specified range (e.g., 65 to 109X PL). This is considered a 2 sigma 
(two standard deviation) response. The shock and vibration loads are used by 
dynamist as input to structural models. 

NASA/MSFC uses similar techniques for developing random vibration criteria for 
the total launch vehicle. A discussion of this approach is found in Reference 
5. 

Chugs are another transient condition in the combustion process that are 
tracked along with the pops. Chugs are low level pressure oscillations whose 
responses on accelerometers would not be discernable. These oscillations are 
obtained from pressure transducer measurements. Chugs will be discussed along 
with the pops since they are both combustion stability type loads and are 
potentially dependent on each other. 

An overall summary of the type of load, information base and typical limiting 
concern is listed in Table 2. 
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TABLE 2 

PROBABILISTIC DYNAMIC LOADS 
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Steady All accelerometers Engine Structural Extensive samples 

State AMS or composite mainstage loads criteria on all 

& PSD plots document accelerometers. 

steady-state 
vibration spec. 


3.4 Mechanical Vibration Loads - General Discussion 


In a rocket engine there are two primary sources of energy which develop 
mechanical and flow vibration loads; these are the combustion process and 
turbomachinery. Rocket engine scaling methodology was developed by Barrett 

and reported in Reference 4. Barrett recognized four sources of excitations: 

1. Mechanical energy from rocket engine fluctuations (i.e. 
combustion) . 

2. Acoustics from the rocket engine. 

3. Aero loads from boundary layer fluctuations. 

4. Self-generating machinery. 

The acoustics and aero loads are primarily vehicle-related excitations, and 
the combustion and machinery loads are more engine-related. 

Barrett's approach to defining scaling parameters is summarized as follows: 

Any structural response possesses a vibration power, P y ., . Likewise, the 
impinging acoustic or flow loads can also be defined in terms of power, 
Pmech’ The two can be related to a vibration efficiency factor, y, as 

P vib 

Y = — 
mech 

This relationship stipulates that a certain portion of the flow, acoustic, 
combustion and machinery power is transferred or absorbed by the component as 
vibrational power. 
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The mechanical vibration power can be expressed as 


P 


vib 


1 

2ir 



cps 


x Wg 


where W = the effective structural weight of the component 

= power spectral density (PSD) of the vibrating 

structure's acceleration 

cps 

g = acceleration due to gravity 

Af = effective bandwidth 

f 

The mechanical power is 

Pmech = TV 


where T is the thrust of the engine and V is the exhaust velocity of the 
rocket engine. 

Substituting in these equations result in 

Y = AfW(G 2 /cps) g 
2irf TV 

Assuming similar structures in different rocket engines possess similar 
dynamic characteristics, the mechanical efficiency factors are equal. 

Therefore: 

(G 2 /cps)n = (jjp) n 

(G 2 /cps)r = (jjp) r 

where r is the reference component, and n is the new component. 
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For the composite or sinusoidal case, the effective bandwidth cancels and 



A mass attenuation factor was also defined where an added mass, W c , is 
mounted to structure where an environment was previously defined. Since the 
mechanical or acoustic forces driving the structure do not change, the 
amplitudes are decreased by a factor of 

H 

W + w c 

The above equation is then modified to: 



TV r 


was assigned a constant value on the component type, and W n + W c was the 
component weight and G r is a PSD of normalized G r vs frequency. PSD's are 
furnished for combustion chambers and turbopumps. 

For example: For the combustion action 



where W n = combustion chamber + nozzle weights 
where W n = fuel + oxidizer weights. 
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The standard approach used by Barrett and still in use today was to define 
the environments by enveloping representative PSD data. Sinusoidal forcing 
functions were only generally addressed by Barrett. 

The Barrett approach is somewhat a broad brush approach in that it only uses 
the gross engine thrust and exhaust gas, TV, as scaling variables for any 
engine forcing function. A more appropriate generic approach is to relate 
the power of each individual energy generating component, e.g. each 

combustor and each turbopump. For instance, the SSME has 7 primary sources 
of energy — the 3 combustors - main injector/chamber/nozzle and two 
preburners, and four turbopumps - HPFTP, HPOTP, LDOP and LPFP. For the 3-2 
engine, there were two combustors - the main injector/chamber/nozzle and the 
GG, ana tne two separate turbopumps, LOX and fuel. For an engine like the 
F-l , the two turbopumps were mounted together with one turbine and 

constitute one turbopump assembly. With the above variations in engine 

components and related engine cycles, it is apparent that generic load 

definitions are best related to classes of components like combustors and 
variations of turbopumps rather than overall engine scaling. The TV 

parameter is inflexible to major engine configuration changes. 

Combustor Loads. The engine main injector/chamber/nozzle environment 
definition can be handled directly by Barrett's method in that essentially 
the total thrust and exhaust velocity are developed by these components. 
The preburners and gas generators, GG's, can be scaled similarly except the 
component injector pressure times area is the T, and the injector velocity 
is the V. Additional sinusoidal loads are superimposed on the random 
levels. These sinusoids are usually turbopump phenomena but can be a 
combustion instability phenomena. 

Turbomachi nerv Loads. The turbomachinery loads are not as easily 
generalized. There are at least two primary load generators: the inertial 

unbalance loads from the turbopump rotor that are primarily related to one 
per rev loading and the flow "noise" loads. The 1/rev is primarily a 
sinusoidal response, whereas the flow noise is primarily a random response 
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plus sinusoidal response at multiples of pump speed. Other sinusoidal 

forcing functions can occur from items like bearing deterioration and 

rubbing. Therefore, the turbopump loads have both a random level and 

sinusoidal components. The sinusoidal components can be transmitted to 
other power generating components like combustors and can be an important 
part of their environment. 

Current turbomachinery sinusoidal load correlation methods essentially use 

2 

Barrett's procedure or pump speed squared, w , as the scaling parameter 

2 

for the composite response, G. The <o can be related to both 1/rev 

2 

rotor loads or flow noise response. So the w scaling can potentially 

be used for both sine and random response. Another approach is to use 

2 

turbopump power (speed times torque) as a scaling parameter for a G 
random response. This would be more in line with Barrett's power and 
efficiency factor concept. The 1/rev and flow noise loads related to flow 
interruption from vanes or impeller blades, etc., are primarily sinusoidal 
responses. 

Vibration Loads - Generic Environment Definition 

The vibration loads for engine components will be defined as both a 
composite and a PSD load function. The PSD will be separated into random 
levels and sinusoids. 

Figure 9 summarizes the approach planned for developing the generic random 
vibration loads. One environment is planned for injector LOX post loading 
using the OPB accelerometer measurements. Another typical environment 
planned is a turbopump environment based on the PBP accelerometer. This 
response is typical of the input response applicable to the HPOTPDD. 
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DEFINED AS PSD ENVIRONMENT 


A 


• SERIES OF SEGMENTS OF G 2 VS FREQUENCY 

• MEAN 

• COV 

• VIBRATION LOADS ARE RESPONSES TO 
TWO BASIC GENERIC LOAD DRIVERS 

• COMBUSTION. I.E. OPB 

• TURBOMACHINERY. I.E. HPFTP 

• FLOW NOISE ■° 01 ' 

• UNBALANCE 

• POWER 



1500 2 TO 20 

khz 


• DEFINE GENERIC ENVIRONMENTS FROM 109% P.L. AND R5 ENVIRONMENT 

• OPB 

• HPFTP 


• USE BARRETT CRITERIA FOR BASELINE GENERIC COMBUSTION SCALING 


• G2 -» F (TV) 
W 


T - THRUST 

W - WEIGHT OF COMBUSTION COMPONENT 
V - EXHAUST VELOCITY 


• MEAN VALUES OF G 2 SCALED 


• FOR TURBOMACHINERY SCALING USE 


• FREQUENCY VARIATION 
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• NO SCALING OF Fj WITH POWER 

. USE 26 TEST PHASE I ENGINE DATABASE FOR MEAN LEVEL 
AND FREQUENCY (G 2 , Fj) BASIS 

• USE R5 UPPER BOUND ENVIRONMENT WITH 26 TEST DATABASE 
MEAN TO 

• DEFINE 2<r& COV'S 

• ADJUST LOW FREQUENCY DC BIAS OF MEAN LEVEL 

Figure 9. Generic Random Vibration Loads 
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The PSD is approximated by a piecewise linear, and its line segments are 
defined based on a mean and coefficient of variance (COV) for both the 
response level, G 2 , and the frequency, F. . The phase 1 engine 26 
test data base will be used for defining (G 2 , F 2 ) points using the 
OPB and PBP response accelerometers. The data has been processed so that a 
mean response level, rather than an envelop, can be defined. A 2a bound 
of the response will be based on the current R5 envelop response for the 
measurements. The available processing of this data has an upper bound of 
2500 Hz. Current plans are to reprocess the basic test data and expand the 
frequency band to 5000 Hz to be compatible with future processing of phase 2 
engine data. More accurate measurements of a 2a or COV response will also 
be available from this reprocessed data. 

The approach for defining the sinusoidal environments is summarized in 
Figure 10. The same 26 test data base and R-5 environment limits will be 
used for their definition. A mean and COV is defined at each discrete 
sinusoidal response. The frequencies are correlated as a function of pump 
speed to allow for frequency shift with power level and for correlation with 
known geometric flow interruptions in the turbopump. The same two engine 
measurements OPB and PBP are used to give data consistent with the random 
response definition. 

The response variation as a function, of power level will be evaluated for 
both types of environments using a power level equivalent of Barrett's 
criteria and as a function of pump speed squares, co , for the 

turbomachinery loads. These two approaches furnish different power scaling 
factors. Figure 11 furnishes an initial comparison of the two approaches 
and shows that w 2 is essentially the power level raised to a 1.3 

exponent, whereas the power function method is essentially linear with power 
level (thrust). The question of transmissibility of the turbopump generated 
sinusoidal responses needs further thought and development. Figure 12 

addresses some of the issues and observations from data. 
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PRIMARILY TURBOMACHINERY RELATED 
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Figure 11. HPFTP Mechanical Vibration vs Power Level Correlation Parameters 
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Figure 12. Sinusoidal Environment 



3.5 SSME Test History Experience and Potential Problems - Pops and Chugs 
The one known significant problem on SSME associated with pops and chugs was 
an ASI line that ruptured during cutoff when the chug - a transient 

pressure oscillation - sucked hot gas and hydrogen into the ASI and ASI 

line. This resulted in a large magnitude detonation or pop. The pressure 

wave from the detonation ruptured the ASI line as noted in Figure 13. The 
fix was to change the shutdown purge operation. 

There are significant variations in pops and chugs test-to-test and 
engine-to-engine. Duty cycle changes like cutoff level and purging are 

variables that affect pops and chugs. During a slow starting engine like 
the SSME, pops occur from both local gas pockets. Occasionally, there is a 
preourner pop into the high power regime (probably late LOX post ignition). 
After cutoff, the pops are the result of combustible gas pockets in the hot 
gas and preburner zones. No known pops have occurred in the main injector. 
The probability of a pop occurring is a function of the start and cutoff 

sequence and length of time during start. On a GG engine system like the 

3-2 that had a faster spin up start, local detonations are probably not 

separable from the basic transient loads. The potential for cutoff pops in 
the enclosed GG hot gas system, though, is there. 

Pops are tracked on the SSME by time of occurrence and maximum peak-to-peak 
magnitude of the pulse. Chugs are tracked by frequency and magnitude. Even 
though there should be some inter-relation between pops and chugs, it is not 
apparent when the tracked parameters are over-plotted, e.g. see Figure 14. 
Reference 5 furnished background information on both pops and chugs from a 
general rocket engine standpoint. 

Pops 

The pop information has been tracked throughout the SSME engine program and 

has recently been translated into a computerized database. Separate files 

have been made for the individual preburners and main injector. The fuel 
preburner data is being used in developing a generic load model and the 
remaining data will be used for validation purposes. 
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SHUTDOWN 



Figure 13. SSME Fuel Preburner 
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Figure 14. SSME FPB Cutoff Pops and Chugs at A-l 


Figure 15 shows example plots of some of this information for the fuel 
preburner (FPB). The first plot relates the peak "pop" magnitude versus 
time and the second plot furnished the number of "pops" versus time. Little 
direct measurements such as high frequency pressure transducer data are 
available to correlate these vibration shock responses to actual engine 
variables, but expert opinion and technical reports will be developed for 
use in generalizing this data. As with many of the other variables on the 
engine, few comparable measurements are available from previous engines. 
The SSME has had much more extensive measurements than earlier production 
engines like the 0-2 or F-l . 


The pop and chug loads used in the development of a typical baseline set of 
i in or Tfiatiori should be available in the expert systems to aid the user in his 
understanding of a specific load model. From a new user standpoint, one 
needs to have: 


1) information defining the load and its cause, 

2) how it is measured, 

3) how it is processed for use, 

4) whether a physical model is available, 

5) potential concern for damage, 

6) type of event, 

7) key variables, 

8) probabilistic model. 


In addition, a baseline set(s) of mean values and coefficients of 
vari ati ons-COVs (or other parameters) are required. 


Table 3 outlines a proposed format for furnishing this information to the 
user of the loads expert system. It covers the points listed above in a 
logical fashion. The baseline mean values and COVs for the load model will 
be added so that a user can judge whether they are adequate for his 
application. 

Most of the information in Table 3 is self-explanatory except for the 
probabilistic model related items. In the probabilistic model, the 
engi ne-to-engi ne variation is considered an independent random variable 
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Table 3 

Pops and Chugs 


1.1 POPS 

1.11 DEFINITION POP A LOCAL DETONATION IN THE COMBUSTION 

ZONE REGION AND/OR ADJACENT HOT GAS 
SYSTEM. POPS OCCUR AT START AND CUTOFF. 

THE APPARENT CAUSE IS EITHER LATE INJECTOR 
ELEMENT IGNITION OR LOCAL POCKETS OF 
STATIFIED GAS THAT REACH A DETONATION CONDJTON- 
APPROP I ATE TEMPERATURE, PRESSURE AND MIXTURE RATIO. 
(STORABLES ALSO HAVE POPS AT MAI NSTAGE * ASSOC I ATED 
WITH INJECTING LIQUIDS AND MIXING) 


1.12 HOW MEASURED 


ACCELEROMETERS LOCATED ON LOCAL STRUCTURE. 
EXTERNAL TO INJECTOR. 


1.13 HOW PROCESSED PEAK MAGNITUDE AND TIME OF OCCURANCE 

WITHIN DEFINED TIME PERIODS. 

SHOCK RESPONSE OF REPRESENTATIVE DATA. 
STATOS RECORDS FOR QUALITATIVE LOOK. 


1.14 PHYSICAL MODEL TBD 


1.15 POTENTIAL CONCERN 
FOR DAMAGE 


START SYSTEM SMALL LINES 
LOX POST/INJECTOR 
TURBINE BLADES AND NOZZLES 
SHEET METAL 

INSTRUMENTATION PROBES 


1.16 TYPE OF EVENT COMMON IN GENERAL * SUFF I CENT SHOCK LEVEL TO CAUSE 

DAMAGE TO SSME , RARE EVENT IN 
STEADY STATE OPERATION ON SSME 

1.2 GLOBAL VARIABLES ENGINE TYPE 

STAGE COMBUSTION- MAIN IN JECTOR , PREBURNERS 
AND DUCTING AND ASSOCIATED MIXTURE RATIO 
-SLOW START BOOTSTRAP ENGINE 

GAS GENTERATOR- MAIN INJECTOR, GAS GENERATOR 
AND DUCTING 

EXPANDER CYCLE- MAIN INJECTOR 
-SLOW START 


1.3 G-VAR I ABLE STAGE COMBUSTION 

1.31 I -VAR PREBURNER(PB) MAIN INJECTOR(MI) 

1.4 1 1 -VAR 

1.41 GENERAL 

1.411 ENGINE TO ENGINE VARIATION I I 

1.412 DUTY CYCLE CHANGES DET DET 

(START & C/O SHAPE & TIME, 

SS LEVEL) 

1.42 START 
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Table 3 (Continued) 
Pops and Chugs 


STAGE COMBUSTION 

PREBURNER(PB) MAIN INJECTOR(MI) 


1.421 CHUG 

1.422 LOX POST IGNITION 

1.423 LOCAL GAS POCKETS 

1.424 8LOWBACK-OPB 


1.43 STEADY STATE 

1 l 7 A I r>v D'”'kC' T I 

1.44 CUTOFF 

1.441 CHUG 

1.442 LOCAL POCKETS MIXED WITH HE I 


1.45 LATE POPS 
1.451 LOCAL GAS POCKETS 


1.5 PROBABILISTIC MODEL 


1.6 DATABASE FORMAT OF LOADS 


START FAC*ENG*DC*(BB+LP+GP) 
SS FAC*ENG*DC*LP 

C/0 FAC*ENG*GP 

POST C/O FAC*ENG*GP 


DISTRIBUTION OF MAXIMUM LOAD IN POP 
cDcrTDiiM ctai pd TO MAXIMUM LEVEL OF POP 


FAC'SCALE FACTOR 

1.411 

1.412 

LP-LOX POST INGtTION 
GP-GAS POCKETS 
I • INDEPENDENT SERIAL 
IP- INDEPENDENT PARALLEL 
DET 'DETERMIN1ST IC 


NOTE: BLOWBACK CAN OCCUR WHEN THE EJECTOR DOME CAVITY HAS SOME GAS IN 
AND BACK PRESSURE PUSHES HOT GAS UP THE INJECTOR TOST. 

IN THE SSME THIS OCCURS IN THE OPB^BUT NOT THE FPB WHICH HAS ITS 
DOME PRIMED TO A LIQUID STATE PRIOR TO THE BACKFLOW STATE. 

J-2 HAD GG BLOWBACK INTO VALVE SEAT, START 
SEQUENCE CHANGED TO ELIMINATE. 


IT 


unT p . THF PROBABILISTIC MODEL CONSIDERS THAT THE TIME PHASING WITHIN 
tSe SPLIT UP PORTION OF THE DUTY CYCLE IS NOT IMPORTANT TO THE COMPONENT 
LOADS E G. THE DAMAGE POTENTIAL IS THE SAME THROUGHOUT THE START TRANSIENT. 
THIS MAYBE CHANGED LATER. 
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that has a wide variation engine-to-engine. Start and cutoff duty cycle 
modifications are more of a deterministic parameter. Pops and chugs are 

different during start, steady state, and cutoff, so their probabilistic 
parameter estimates are also different. Preburners and main injectors also 
have differences, so they are also separated. Pops and chugs typically 
occur during transient conditions, but rarely occur during steady state 
operation. 

Chug Loads . The chug load format has been prepared in keeping with the pop 
data approach for presenting data (see Table 4). 

The text information is self-explanatory. In this case, a physical model 
can readily be developed for use in the options portion of this, contract. 
Chugs occur each test during start and cutoff, and like pops are a form of 
combustion related instability. The probabilistic model is simpler than 
that proposed for pops since there is not the randomness of occurrence nor 
the multiple causative variables. 

I nternal Flow Dynamic Loads . Internal flow (i.e. inside ducts or components 
rather than external flow like shell flutter) has become critical 
environment on high energy flow systems like the SSME. A series of problems 
has occurred throughout the engine development that are related to fluid 
dynamics. The major problems have been reported elsewhere in the literature 
by Rocketdyne and NASA. A good summary of both vehicle and rocket engine 
related problems are summarized in Ref. 6. From Table 5 (reproduced from 
that document), it is readily observed that environmental problems on 
earlier engines were not a problem. 

The only one noted was in 3-2 (APOLLO) engine where engine propellant line 
bellows lack fluid structure coupled vibrations and failures. There are 8 
listed for the SSME engine. These problems have spawned extensive analysis 
and testing within Rocketdyne and NASA to understand the associated 
phenomena and started toward developing better predictive environment 

methodology. All four components in this project have significant fluid 
dynamic loads. 
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Table 4 


5.5 


5.51 CHUG 

5.52 DEFINITION 


CHUG 


SSME 

ci:u: 


S1NUSIOOAU WAVE SHAPE THAT GROWS IN A LINE OSCILLATION IN ONE 

« ffi®>™rS"SE‘SK SSSSM TSIns.ent 

™c“u “ ml 'll SSSIy&Y Kll UNDERSTOOD, STEADY STATE C«U0 
99 N ?J.I*,°.;!c hirciv THAN TRANSIENT CHUGGING. 




5.53 HOW MEASURED 

5.54 HOW PROCESSED 

5.55 PHYSICAL MOOEL 

5.56 POTENTIAL CONCERNS 
FOR DAMAGE 


5.57 TYPE OF EVENT 


5.6 GLOBAL VARIABLES 


PRESSURE TRANSDUCERS LOCATED CHAMBER CAVITIES- PREBURNERS, 
GAS GENERATORS OR MAIN CHAMBERS. 

1 SPJ 1 asffssra.swnas - T “ — 




rHIS FORM OF OSCILLATION J HAY T DO [|0 “AMAGE AT ALL^ ^ joints 

S'S fSlo**"™ 1 MdSedS'oF 'oSCILLAT.ONS Of VARYING 
uifitiiTiince DCO TRANSIENT. 


(ARE FOR STEADY STATE CHUGS. 

*«!«£ KJS'SS.RSS-.y start and CUTOFF. 


ropeuaht feed SYSTEM .motor D|LTA pressure fluid inertanoe, 

I«DIIY 


5.7 PROBABILISTIC VAR IAIBLES 


5. 

5. 

5. 


■NGINE TO ENGINE VARIATION 
1UTY CYCLE CHANGES 
>Ut$E PARAMETERS- 
4AGN ITUOE f DURATION ( 
FREQUENCY AND 
4UMBER OF BLOSSOMS 


I 

DET 

I 


•SCALE FACTOR 
•5.71 

SM^NUMBER OF BLOSSOMS 
JR-TIME DURATION 
:Q-OSCILLATION FREQUENCY 
'.u.ppav: MAGNITUDE OF BLOSSOM 


5.8 PROBABILISTIC MOOEL 


START 

SS 

CUTOFF 


FAC*ENG*DC*<NBLSM,TDUR,FREQ,MAGN> 

II 
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Table 5 

Dynamic Problems Experienced in Flight Vehicles 



The HPOTPDD high pressure oxidizer turbopump discharge duct, the fourt 
component in the CLS study, is in the oxidizer pump discharge system where 
the flow environment is a key load, see Figure 16. High frequency pressure 
measurements such as those shown in Figure 17 are available. Also 
considerable study is currently under way related to a problem the 
injector LOX inlet area of the SSME. Table 6 lists the oxidizer system 
fluid flow problems as well as other hot gas system problems and Figure 18 
shows three of the key variables that influence the loads in the oxidize 
system The power to weight ratio, pump pressures and dynamic velocity hea 
Al a'll doubled or tripled relative to other flight engines. The single 
variable that most effects the fluid-structural Infraction In he her wa 
is probably the high value of the pump velocity head. In the SSME HPOTPDD, 
his LL, is much greater than on any other Rocketdyne tor opump or 
engine. Generic flow loads studies as part of the duct loads should aid 
setting better limits of flow related parameters in new hardware designs. 

Historically, conceptual or preliminary sizing of engine components are done 
scaling previous engine components using strength parameters, not 
environmental loadings like vibration or flow. The as effort shoud 
furnish additional criteria to make a more accurate sizing assessment 
starting from a conceptual design standpoint. 

Similarly, the hot gas system has had a series of problems where the other 
three components under study are located. The initial fluid flow loading 
being addressed in the hot gas system components to support the transfer 
duct and LOX post load modeling. Host of the SSME fluid 
modeling and measurements have been done in this system. The HPOTPDD 

modeling will then be developed. 
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SSME POWER HEAD 















TABLE 6 

SSME FLOW AND FLUID STRUCTURAL PROBLEMS 


Problems in HPOTP Discharge System 

• Flow Straightener (for Flowmeter) 

• Flow Meter 

• Oxidizer Valve 

• HPOTPDD Lip Cracking 

• Main Injector Inlet Vane - 4000 Hz 


Problems in Hot Gas System 

• Hot Gas Manifold Flow - Fuel Side 

• Fuel Transfer Duct Coolant Liner - Fuel Side 

• Preburner Lox Post - Fuel Side 

• Main Injector Lox Post 

• Bellows Shield - Fuel Side 

• Turbine Blades 

• Nozzle Steerhorn - Start Transient 

• ASI Orifice - Cutoff Chug and Pop 

• HPFTP Kaiser Hat/Nut Failure 

• Temperature Probes 
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Figure 18. Key Factors for High Pressure Oxidizer 
System Flow Environment 



4.0 TECHNICAL PROGRESS AND PROBABILISTIC MODEL DESCRIPTIONS 
4.1 Introduction 

The development of the probabilistic model for composite load descriptions 
during the second year has focused on the following topics: 

(1) Data base development 

(2) Transient load modeling 

(3) Mission phase modeling 

(4) Improvements to the probabilistic model 

(5) Inclusion of expert opinion data 

(6) Generic load calculations 

The data base development has been progressing throughout the program. This 
is one of the essential developmental areas, since it is where the primary 
interaction between the probabilistic model and the expert system takes 
place. 

The transient model was developed because the treatment of the loads during 
a transient event is fundamentally different from the treatment of the loads 
during the quasi-steady and steady state portions of the mission. In the 
transient model, the peak amplitude and the time of occurrence of this peak 
are treated as random variables. The nominal behavior is not separated from 
the load variation in this model. 

Because there are up to 42 individual load and engine parameters which can 
influence the composite load calculation, it is important to provide a 
continuous, realistic transition between the three mission phase types. 
This linking of the mission phases has been accomplished, and compares well 
wi th avai labl e data. 

An improved version of the primary probabilistic model, RASCAL, was 
incorporated into the computer code system. This version of the program 
allows the user to direct importance sampling schemes via input. Therefore, 
loads which are rarely occurring, but are potentially important for design 
or failure analysis, can be examined quickly. 
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Data from SSME designers and analysts on their expert opinion about the load 
variability was obtained and incorporated in the probabilistic model and 
data base. Additional information will be added as it becomes available. 

Late in the year generic engine calculations were performed. These 
calculations will be improved and updated, as well as proceeding with the 
validation and verification, during the third year. 

The following sections provide additional details about each of these areas. 
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5.0 PROBABILISTIC LOAD ANALYSIS FOR GENERIC SPACE PROPULSION ENGINES 

5.1 Introduction 

The development of a probabilistic load model for a generic space propulsion 
engine has been proceeding in several steps. At this time it is wise to 
reexamine these steps to illustrate how the program encompasses the goal of 
coupling the probabilistic load model, which predicts how mission and design 
changes will affect the critical loads in the specified engine, with the 
development of an expert system for load prediction. 

5.2 Probabilistic Models For Generic Engines 

To examine in detail how the probabilistic model will deal with generic 
engines, it is necessary to examine first the relationship between the 
probabilistic model and the expert system. Following this study, the 
developments in specific parts of the model during FY86 will be presented. 
Next, examples of the use of the model and the validation work to date will 
be presented. Appendix A furnishes details about using the code in a stand 
alone mode used during development. The probabilistic code ANLOAD (ANalyze 
LOADs) is being incorporated into the expert code by Rocketdyne. 

An overall picture of the flow of information is provided in Figure 19. 
This is not meant to represent the current status of the expert system being 
developed by Rocketdyne, but rather is a representation of the information 
flow between the probabilistic model and the expert system. Some detailed 
discussion of this figure is warranted. 

The critical information which must be communicated between the 
probabilistic model and the expert system is the mean, variance and 
distribution type for the individual loads and for the composite loads. As 
additional analyses are ' performed the database of information will be 
updated and previous analyses will be saved for future requests. Therefore, 
the important development work is in generating the new probabilistic 
information for individual and composite load parameters. 
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Figure 19. Interaction Between Probabilistic Information and Expert System 
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Figure 19. (Continued) 
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There are three boxes identified as tables; specifically, tables for the mean 
and variance coefficients, and the type of distribution which describes the 
data. In actual practice, these are not look-up type tables; rather they are 
function tables, e.g., those for the SSME influence functions. For each of 
these tables, the entries can be based on either previous data or analysis, or 
on scaling of data representative of one engine type to another. For example, 
the turbine speed in the 02 engine may not be known but, given the power 
requirements, it may be reasonably approximated by scaling known SSME data 
according to the power requirements. Thus, when individual parameters are 
unknown they are estimated by providing scale factors from better known, or 
understood, data bases. 


These scale parameters must also be provided for both the mean and variance 
coefficients. One can envision a situation in which the mean value will scale 
based on one or more parameters, while the variance will be scaled based on a 
different set of parameters. For example, returning to the example of 
calculating 02 loads based on SSME loads, the turbine torque may scale 
according to the horsepower and speed ratios, but the variance will have to 
account for other differences in the engine. The reason is that the variance 
in the turbine torque is dependent on the other such variables. For example, 
head rise split between the two sequential pumps, the basic engine control 
philosophy, and the differences in a gas generator driven and a stage 
combustion cycle. Therefore, the table of variances must not only contain the 
total variance for the individual and composite loads, but it also must 
provide information on how this total variance is partitioned among the 
individual engine parameters and/or loads. 


Once such tables have been defined, the expert system can then begin to quiz 
the user on the engine type and mission profile to be examined in the current 
analysis. This is identified as the Composite Load Request. The expert 
system then decides which are the critical loads for the problem to be 
analyzed and selects from the table the appropriate mean and variance 
coefficients, together with their associated scale factors. If a request for 
an analysis has been made in which all of the individual loads are known from 
previous analysis, then the scale factors are all set to one. If new 
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individual load data must be generated, the code selects that engine which 
is most closely aligned with the requested engine type. A specific turbine 
is also selected to reflect differences in the oxygen and fuel sides, and 
the appropriate scale factors are chosen. 

Assuming that this is not a composite load analysis that has been performed 
previously, each of the individual loads to be included in the analysis are 
checked to insure that they are characterized probabilistically by their 
mean, variance and distribution type values. If the distribution type is 
non-normal, then the appropriate transformation is selected to calculate the 
distribution parameters based on the type of distribution which describes 

the individual load. Currently, it is assumed that the coefficient of 
variation for the individual load, as well as the composite load, is 
independent. If this is not the case, a quasi-steady analysis is called for 
in which the distribution parameters are allowed to vary in a time dependent 
fashion. 

At this point all of the necessary probabilistic information has been 
collected and the probabilistic synthesis of the data can be performed. 
This synthesis is done using one of the three probabilistic models: (1) 
Monte Carlo, (2) RASCAL, or (3) QI_M. The results are then sent to a post 

processor for display. Finally, the results of this analysis are placed in 
the data base for future reference. 

5 . 3 Link i n g Pi f fere n_t_M i ssion Histor y Phases 

The insertion of a "probabilistic" model in the flowchart of Figure 1 is an 
over simplification of the actual process of interaction taking place 
between the expert system and the probabilistic analysis. In the actual 
analysis the expert system must be constantly updating the input to the 
probabilistic model so that the appropriate techniques, and data bases are 

used. One of the more difficult transition regions in which to perform this 

link occurs during engine start-up and subsequent power up to the demanded 
thrust levels. In this type of analysis three different mission 
requirements are demanded of the probabilistic model: (1) transient 
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analysis, (2) quasi-steady analysis, and (3) steady state analysis. 
However, each cannot be performed independently of the other since the loads 
are continuous functions. Thus, it would be inappropriate to have the 
transient analysis predict the (mean) temperature at the end of the 
transient phase to be 1000°R, while the subsequent quasi-steady state 
analysis is predicting a temperature at the start of the quasi-steady state 
analysis (i.e the end of the transient phase) to 2000°R. Thus, the 
transient analysis must be able to predict the load behavior during the 
defined time period, as well as provide a smooth transition to the 
subsequent mission phases. Similar arguments apply to the quasi-steady and 
steady state analysis. However, the difficulty in these situations is eased 
greatly when there are adequate functional relationships between the 
different phases, for example, as in the case for the influence functions 
for the SSME engine. 

To illustrate the method for dealing with the transient response, an example 
using the SSME HPFTP temperature was examined. Figure 20 shows the analyses 
of three engine tests performed by Rocketdyne to calculate the turbine inlet 
temperature based on a combination of engine measurements and a turbopump 
model. Tests 902349 and 902363 show three distinct peaks (The third peak is 
much smaller in magnitude than the other two and occurs near 1.8 seconds.) 
in the temperature, while test 902356 appears to have only two peak values, 
at least relative to the other two tests. In addition, these peaks occur 
over a relatively narrow time period, on the order of tenths of seconds. 
The two questions to be addressed are: (1) how should the variable number 
of peaks be handled?, and (2) how should the variable magnitud e of the peaks 

be handled? 
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Transient Temeperature Model 



Figure 20 Selected Tests For the Temperature In the SSME HPFTP 
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5.4 Steady State 


The input to the AN LOAD program requires that the beginning time and power 
level as well as the end time and power level be specified by the user. In 
some cases one or more of these parameters is forced by the program to 
ensure a continuous time line and power profile. During the quasi-steady 
and steady state phases of the time history, the influence coefficients are 
used to determine the loads seen at the critical components of the engine. 
During the transient phase, only the peak load value and the mean time of 
occurrence of this peak, amplitude are of concern. 

The testing of this transient model has been performed. There were several 
problems with the implementation of this method, the most significant one 
being the assignment of the variability of the independent engine parameters 
below the 65% power level. Because several of the parameters have 
non-physical predictions for their values below 65t power it is necessary to 
restrict the quasi-steady state phases to be applicable only above this 

power 1 evel . 

5 . 5 Transient Load Model Deve lopment 

The computer program PEAKS has been constructed to create a response 
envelope from the observed transient responses of the input variables to the 
influence functions. The response envelope defines the beginning, apex, and 
end of the individual transient events both in the magnitude and time 
domains. These critical points, the start, apex, and end, are hereafter 
referred to as knots. Due to the short duration of a transient response, 
random variation between consecutive knots is neglected, i.e. the response 
is assumed to be piecewise linear between knots. 

PEAKS has the option for selecting various levels of accuracy when it 
constructs the response envelope. The first level of accuracy is provided 
by examining the first derivative of the response function to determine when 
the peak values occur. The calculation of the derivative can be made using 
the standard finite difference approximation, using central differences, or 
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can employ a four to twenty point moving average. This is done for the 
situation in which the data are oscillating about a mean trend line that is 
monotonic. In such a situation, the correct determination of the knots 

requires that these oscillations be smoothed to some extent so that true 
peak values can be observed. Example plots of these predictions are given 
in Figures 21 through 23. 

PEAKS also has the ability to examine the second derivative (with the 
central difference finite difference approximation) to further refine the 
selection of the knot points. 

The transient model for the example being considered is shown in Figure 23. 
In this figure the "+" symbols represent the average response of the data. 
The open boxes represent the mean transient model response. Finally, the 
solid lines are the bounds representing the two standard deviation spread 
about the mean value. As can readily be seen from Figure 21, if the lower 
bound is used, there will be only two peaks in the analysis. 


Transient Temeperature Model 



Time Afler (teeo'wdt) 

C TrcnsleM Mcdet * Average 


Figure 21 Probabilistic Model For Transient Analysis 
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The linking of the transient model with the quasi-steady state model should 
be checked at this point. For the SSHE. the data analysis for temperature 

ends at approximately 5.0 seconds. From the data analysis, this appears to 

be at approximately the 86X power level. Averaging the available test data 
and calculating the standard deviation gives a mean and bounds for the 

temperature response. Then, this is compared with the calculations obtained 

from the influence function in order to calculate the temperature. The,e 
results are shown in Figure 22 where the break in the plot represents the 
break between the data analysis (less than 5 seconds) and the predictions 
using the influence functions (greater than 5 seconds). 


Transient Temeperature Model 



Figure 22 Transient And Quasi-steady Model Interaction 
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As this plot clearly shows, the link between the actual, observed data for 
temperature and the predictions between the temperature by the influence 
functions is quite good. To illustrate this point. Figure 23 shows the mean 
and bounds for the temperature if the transient data had simply been 

extrapolated from the 86% power level up to the full power level of 104%. 

The thick lines represent the results of the extrapolation. Clearly, the 

transient model is describing the temperature behavior well. Thus, the 
transient model described here provides the appropriate model for the 

transient and transition from transient to quasi-steady state analysis. 

These models, and previous analyses with the steady state analysis, have 

shown that reasonable, cost effective, and accurate results can be obtained 
for space propulsion engines. However, because of the current information 
in the data base, almost all of the verification and validation of the 
computer codes have been performed for the SSME. There still must exist the 
capability to address, not only current engine, but also the model must be 
able to account for mission operations outside of the current experience as 
well as design changes. Clearly, radical departures from the existing 
engine types will be less accurately handled by the expert system, however, 
a capability for perturbations on the present state of knowledge should be 
manageable in the computer model. 

5.6 Database Development 

The development of a standard database for generic space propulsion engines 
is included in the computer code system. The current version of the code 
incorporates the data which have been derived from analysis of the 
independent load parameters. The program has been designed to provide the 
user with default values if none are provided during the interactive input 
session. Each time a default parameter is obtained from the data base it is 
identified, as well as the source of the data value from which it was 
obtained (e.g. SSME data analysis, expert opinion, etc.) Each of these 
default values, if selected, is automatically inserted into the input file 
being constructed by the user. 
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The database for generic space propulsion engines has been tested with a 
standard input problem to assess its use for a duty cycle calculation. A 
problem has been constructed which ramps up to 104% power level from 65%, 
remains at this steady state condition for several seconds, throttles down 
to 65% power level, operates at 65% for several seconds, and throttles back 
to 104% power. The current version of the code approximates the data which 
has been used in standard mission history profiles. 

5.7 Using The Current D ata Base 

There are two situations to consider when developing a data base for use in 
a generic type of analysis. First, the analysis may ask for current engine 
designs, or mission requirements currently within the design specifications 
of a specific engine type. For example, one may wish to examine the effect 
of operating the SSME fuel turbopump at 106% power instead of 104% power. 
In these types of analyses, the deterministic analysis covers the range of 


Transient Tenieperature Model 



Time Mler Engine (tecondt) 
* A*ercge Bc-.tic* 


Figure 23 Extrapolation Of Transient Model To Full Power Level 
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physical parameters and the engine performance data and associated 

probabilistic information can be interpolated to predict the range of loads 
in the engine. Currently, the default for interpolation of the 

probabilistic data is to assume that the coefficient of variation (defined 
as the standard deviation divided by the mean) is constant and that the 
shape, i.e. distributional form, of the random behavior is constant. In 

this case the underlying deterministic model is used to obtain the new 

nominal level and the probabilistic parameters are readjusted to agree with 
the assumed value for the coefficient of variation. The probabilistic 
modeling then proceeds as described previously. 

A more difficult use of the data base involves analyses for which no data 
for thai specific engine or mission profile exists. Yet .to be a truly 
generic analysis capability such situations must be addressed. In 
extrapolating the current experience with space propulsion engines to other 
mission requirements or design modifications it is necessary to make some 
assumptions about what will remain constant and what must be changed. The 
following hypotheses will form the foundation for the extrapolation 

process. These are default assumptions. If a user wishes to change one, or 
more, the capability will be provided. 

H ypothesis I. The variable to be extrapolated will be described by the 

same distributional form describing the load variable in the current data 

base which is closest, in a physical modeling sense, to the variable to be 

extrapolated. 

H ypothesis II. The variable to be extrapolated will be described by the a 
scaled value of the mean value for the load variable in the current data 

base which is closest, in a physical modeling sense, to the variable to be 
extrapolated. The mean value will be allowed to vary with the power level. 

H ypothesis III. The variable to be extrapolated will be described by the 

same coefficient of variation (COV) for the load variable in the current 
data base which is closest, in a physical modeling sense, to the variable to 

be extrapolated. The COV will be allowed to vary with the power level. 
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H ypothesis IV. Non-normally distributed random variables will obtain a 
mean value from the appropriate scaling parameter under Hypothesis II and a 
variance from the COV under Hypothesis III. These parameters will then be 
transformed according to the distributional form given by Hypothesis I. 

The manner in which these hypotheses are used will become clearer in the 
following section showing their use in an example. For now, an examination 
of the current data base and the behavior of the variance coefficients is 

discussed. 

Moa n Coefficients- The determination of the magnitude of the mean value of 
variables not in the data base will be handled primarily by the expert 
system. For example, if it is desired to predict the torque for an engine, 
a probabilistic model has no means for estimating what the nominal value 
should be if there are no data available for analysis. If the case under 
study is adding data to the data base, then the data analysis section of the 
probabilistic model will perform the appropriate calculations to estimate 
the mean value. Otherwise, it is assumed that the expert system provides 
the nominal levels for the necessary variables. 

Variance Qpf f i ci ents . While the mean levels of engine system related loads 
are primarily a deterministic quantity, the variability about the mean level 
is primarily a probabilistic quantity. However, unlike the mean, the expert 
system code may intercede to change the calculations described here based on 
data contained in its knowledge base. Thus, it must be remembered 
throughout this discussion that the rules applied here are generic in nature 
and may be modified as deemed necessary by the expert code. 


There are two cases to consider in defining the variance coefficient table; 
(1) the individual load and/or engine parameter contribution to the overall 
variability in either the individual load or composite load being 
constructed, and, (2) the time dependent behavior of this variability. The 
first topic is important because of the differences which can exist between 
engine designs. For example, in the SSME significant variability exists in 
the high pressure inlet pressure. But since it is a two stage engine, such 
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variability is partially caused by the variability in the low-pressure 
outlet pressure. In a one-stage engine, one could imagine significantly 
less variability in the inlet pressure simply because there is no 
intermediate turbine. If a table is available which partitions the 

variability in the load variable of interest among all if the input 
parameters and loads in the calculation of this overall variability, the 
individual components can be "turned" on and off as demanded by the expert 
system. For the current data base, the engine system related loads are 
given by the influence functions: 


L„- I 


i=l 


V L i 


( 1 ) 


L. = Y a. . 

l ^ i,l 

k=l 


Var(L p ) = 


PL 


k-1 


(C0V(L i ))‘ 


Var(L-) 


a i V a r ( L i ) 
M 


■ z b, I (c, , ..pl'- 1 ) 2 . 

i j = l J [ k=l J,1 ’ k 


(COV^X.) 


(2a) 


(2b) 


L.j - nominal engine value for dependent load 

a i,k = nominal engine coefficients for calculating mean 
value of L- at PL. 

PL = power level 
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r = influence coefficient set 

j i 5 k 

COV(xj) - Coeff. Of variation of independent load j 
C0VtL 1 )- Coeff. of variation of dependent load L i 
b ’ , On-off function for the independent variable j 


Equation (5) provides a means for partitioning the variability in the load 
among the independent variables Xj and L, , assuming that the independent 

variables are independent. This table has been constructed and is used in 

the examples discussed below. 

t i mo npnpndent Behavior It is assumed that there is a one-to-one 

correspondence between time and the power level; therefore the subsequent 

discussion will talk of the COV as a function of the power level as opposed 

to time. 


Figure 24 and 25 display plots of the COV for selected SSME variables as a 
function of power level. The symbols represent the calculated COV at one 
percent power level increments, while the solid line represents the best fit 
curve from a linear, exponential, logarithmic, and power curve functional 
forms. Each of the 20 turbine load variables was described with a curve 

similar to that shown in Table 7. In many cases the regression coefficient, 
r 2 was very close to 1.0 and, therefore, a very good fit was obtained. 

However, in some cases, a low value of r 2 was obtained. An example is the 
HPFTP torque in which r 2 is equal to 0.3253. The best fit line together 

with the influence function predictions are shown in Figure 26. Obviously, 
the fit is very poor. The important point to note is that the absolute 
magnitude of the COV changes only in the third decimal place over the entire 
applicable range of the influence functions. In this case, the expert 

system has the model assume that the COV is constant. Based on the data 

analysis, Table 8 presents the nominal levels for the COV as a function of 

the power 1 evel . 
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coefficient or variation 


4 .&% 


SSME HPOTP Inlet Temperature 



Influence Function 


POWER 

+ Logarithmic FH 


Figure 24. SSME HPOTP Inlet Temperature COV as a Function of Power Level 


The amount of information and how it all fits together is best illustrated 
by way of examples. Therefore the following section discusses some 
calculations which have been constructed to illustrate the methods discussed 

previously. 
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COEFFICIENT OF VARIATION 


SSME LPOTP Discharge Pressure 



Figure 25. SSME LPOTP Discharge Pressure COV as a Function of Power Level 


5 . 8 Improvements To The Proba b ilistic Modeling Code 

The predictions currently being made by the AN LOAD program are based on 
simulation methods which are time consuming to produce. The primary reason 
is that low probability events are of interest in this program and the Monte 

Carlo method is slow for such predictions while the barrier crossing method 
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TABLE 7 

COEFFICIENT OF VARIATION AS A FUNCTION OF POWER LEVEL 



LPOTP Torque 

Exponential 

LPFTP Torque 

Logarithmic 

HPOTP Torque 

Power Curve 

HPFTP Torque 

Linear 

1 r- l _ j _ 

Lru i r r I On i q lc 

Power Curve 

LPFTP Flowrate 

Power Curve 

HPOTP Flowrate 

Logari thmi c 

HPFTP Flowrate 

Logari thmi c 

LPOTP In Press 

Power Curve 

LPFTP In Press 

Exponential 

HPOTP In Press 

Linear 

HPFTP In Press 

Logari thmi c 

LPOTP In Temp 

Linear 

LPFTP In Temp 

Li near 

HPOTP In Temp 

Logari thmi c 

HPFTP In Temp 

Exponential 

LPOTP Out Press 

Power Curve 

LPFTP Out Press 

Exponential 

HPOTP Out Press 

Logari thmi c 

HPFTP Out Press 

Logari thmi c 

Linear Curve: 

C0V(P L ) = 

Exponential Curve: 

C0V(P L ) = 

Loaari thmi c Curve: 

C0V(P L ) = 

Power Curve: 

C0V(P L ) = 


-5.212E-01 

1.030E-02 

9 . 999E-01 

1 . 322E-03 

1 . 250E-02 

9 . 985E-01 

-4 . 362E-01 

5.798E-03 

9 . 646E-01 

2.450E-05 

4.675E-03 

3.253E-01 

-C . 293E-01 

3. 170E-03 

9 . 634E-01 

7 . 746E-03 

9.409E-03 

3.334E-02 

2 . 435E-03 

1 .398E-02 

9.941 E— 0 1 

1 .125E-03 

1 . 187E-02 

9.919E-01 

6.598E-01 

3.844E-03 

1 .000E+00 

8 . 889E-02 

4.846E-03 

8.084E-01 

3.366E-03 

3.852E-03 

9. 995E-01 

1.981 E-03 

6.537E-03 

9 . 728E-01 

2 . 295E-04 

7.172E-03 

9 . 964E-01 

-1 .318E-03 

8.632E-03 

8.080E-01 

1 .618E-02 

4.321E-02 

9.986E-01 

-1 .246E-01 

2.701 E-02 

7.072E-01 

-5.095E-01 

4.41 9E-02 

9. 970E-01 

5 . 398E-01 

9.023E-04 

9 . 1 32 E— 01 

8 . 1 28E-04 

1 . 863E-03 

9.991 E— 0 1 

1 .223E-03 

3.285E-03 

1 . 000E+00 


P L + b 
e ax , x = P 
+ a- 1n(P L ) 

n 
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C^F^CI^T £F c 


SSME HPFTP Torque 



POWER 

Influence Function + Uneor Fit 


Figure 26 SSME HPFTP Torque COV as a Function of Power Level 




TABLE 8 

COEFFICIENT OF VARIATION AS A FUNCTION OF POWER LEVEL 


Variable: SSME 

Curve Type 

a 

b 

r 2 

LPOTP Torque 

Exponential 

-5.212E-01 

1 .030E-02 

9 . 999E-01 

LPFTP Torque 

Logari thmi c 

1 .322E-03 

1 .250E-02 

9. 985E-01 

HPOTP Torque 

Power Curve 

-4 . 362E-01 

5 . 798E-03 

9 . 646E-01 

HPFTP Torque 

Constant 


4.701E-03 


LPOTP Flowrate 

Power Curve 

-6.293E-01 

3. 1 70E-03 

9.634E-G! 

LPFTP Flowrate 

Constant 


9.451 E-03 


HPOTP Flowrate 

Logari thmi c 

2.435E-03 

1 .398E-02 

9.941 E-01 

HPFTP Flowrate 

Constant 


1 .189E-02 


LPOTP In Press 

Power Curve 

6. 598E-01 

3.844E-03 

1 . 000E+00 

LPFTP In Press 

Constant 


5.281 E-03 


HPOTP In Press 

Li near 

3 . 366E-03 

3.852E-03 

9 . 995E-01 

HPFTP In Press 

Constant 


6.563E-03 


LPOTP In Temp 

Constant 


7.401 E-03 


LPFTP In Temp 

Constant 


7.367E-03 


HPOTP In Temp 

Logari thmi c 

1 . 61 8E-02 

4.321 E-02 

9 . 986E-01 

HPFTP In Temp 

Constant 


2.400E-02 


LPOTP Out Press 

Power Curve 

-5 . 095E-01 

4.41 9E-02 

9. 970E-01 

LPFTP Out Press 

Constant 


1.533E-03 


HPOTP Out Press 

Constant 


1 .863E-03 


HPFTP Out Press 

Constant 


3.287E-03 


Linear Curve: 

C0V(P L ) = 

a * P|_ + b 



Exponential Curve: 

C0V(P L ) = 

b-e ax , x = P|_ 



Loaarithmic Curve: 

C0V(P L ) = 

b + a • 1 n C P L ) 



Power Curve: 

C0V(P L ) = 

b-PL 
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and QLM technique, unless all of the inputs are normally distributed, are 
too approximate in nature to provide reasonable results. To examine the low 
probability events, either an importance sampling scheme for the Monte Carlo 
method or an improved version of RASCAL must be used. Of course a fourth 
probabilistic method could be employed. This has not been done for reasons 
which were given in the literature review completed earlier. The RASCAL 
method offers a variety of advantages for use in the expert system code - 
one of the primary ones being the ability to have the user define the range 
of the input probability density functions which he wishes to use. This 
capability allows the user to specify, by input, an importance sampling 
method. This capability has been included in the current version of AN LOAD 
and requires no new inputs, since the three parameters needed for input 
currently have been modified so that the third parameter specifies the lower 
limit of the input range to be examined. For example, if the following 
input is made for the fuel inlet total pressure. 


ID Mean Standard Deviation 

2 28.5545 7.38417 

a normal distribution (ID = 2) with a mean value of approximately 28.6 and 
standard deviation of 7.4 is used to describe the random variation in the 
fuel inlet pressure. The value of P 3 set to 0.001 implies that the ^nput 
probability density function will cover the range from the 0.1 to 
99.9 th percentile values. This is opposed to the equal probability 
version which would cover only the range from 100C1/N) . to 
100(1-1 /N) th , where N is the number of bins used to describe the input 
variables DPD. In the sample runs constructed in the calculations shown in 
Figures 27 and 28 presented later in this section, N was set to 20, thus 
the possible range of input values is from the 5th to 95th percentile. As 
the figures show the variability of the results is changed significantly by 
limiting the range of the input distributions. It must also be noted that 
the values shown in these figures are for illustration only, since the 
physically realistic range of the input values was not limited for these 

cal cul ati ons . 




0.001 
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Additionally, a new random number generator was incorporated into the 
probabilistic code which should increase the period of the pseudo-random 

number generator by more than an order of magnitude (to approximately 
17 , 000 , 000 ). 

5 . 9 Qu ick Lo ok Model 

In some analysis only an approximation to the variability of the load is 
needed. In such a case the relatively long running time of the RASCAL or 
Monte Carlo simulation models is not justified. To provide a program which 

quickly calculates such an approximation the Quick Look Model (QLM) was 
developed. 

The basic assumption made in the QLM model is that all of the individual 
loads and engine parameters used to predict the individual and composite 
loads are normally distributed. In this case the influence function tables 
can be used directly to calculate the mean and variance of the output. If 
there are dependencies among the variables then some modification to the 
current program is needed. However, if the correlation coefficient is 
provided, or calculated, then exact solutions are still available. The 
basic formulas used to perform these calculations are given by the algebra 
of normal distributions presented below. In these formulas p represents 
the mean, or expected value of the random variable, and a is the standard 
deviation, i*e. the square root of the variance. 


These formulas are used in conjunction with the influence equations to 
provide the mean and variance estimates of the load variables. Since the 
influence functions currently in the probabilistic load model do not involve 
any divisions all of the formulations are exact (assuming independence), if 
the probability density functions are all gaussian. 
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Statistic Of The Sum: Z = X ±_Y 


E[z] = Ht + h 
2 2 

V[z] = <7 X + o y + 2/?a x a y 


( 3 ) 


Statistics Of The Differencej — Z = X - Y 


E[z] = Mx - h 
V [z] = a, + cT y - 2 />a„a y 


( 4 ) 


Statistics Of The Produ ct Z = X x ,_Y 

E [Z] = \L X x H y + /><7 x <7 y 

22 22 2 2 222 
V[z] = H x Cy + < 7 ,a y 2p(l x [lyO x O y + p o x o y 


( 5 ) 


Statistics Of The Quotient Z = Y /_X 


h_ 1 

<7« 

1 + — 1 

£ - ^ 

/ ** 
i A + 3 — 

Mx 

<7x1 

^x 

J 

1 1 

v[z] 

* jj r 

a x + a y 

ip^fi 4 

Mx 

Mx 2 

n 

MxMy 


<V 

Mx 4 


( 6 ) 
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Two options exist in the computer code for using the QLM model. If the user 
requests that the QLM model be used and all of the input distributions are not 
normal, then the corresponding mean and variance are calculated by the 
appropriate moment transformation. On the other hand, if the user does not 
request the QLM model, yet all of the input distributions are gaussian, then 
the QLM model is substituted. The QLM substitution is made since there is no 
reason to run a simulation to approximate an answer which can be obtained 
exactly with the QLM model. Figure 27 and 28 show comparisons of the QLM 
model to theory and the simulation methods. 

5 -l° Examina tion Of Model Suitability For Low Probability Calculations 

Tne use of the probabilistic load model for the prediction of low probability 
events, such as pops (small, localized explosions in the engine), raised the 
question: are the available probabilistic methods the most suitable for 

addressing these types of calculations? To answer this question a study was 
performed which examined the use of a fast probability integrator, 
specifically the Chen-Lind (C-L) algorithm, as programmed by Wirsching and 
Wu and the RASCAL program. 

The C-L algorithm is an extension of the technique originally proposed by 
Rackwitz and Fiessler, Ref. 8. In this methodology, the Hasofer-Lind safety 
index is the value of the response variable that minimizes the distance from 
the failure surface to the origin in normal probability space, and the slopes 
of the probability density function are equal. The probability of failure is 
approximated as the value of the normal cumulative distribution function at 
the Hasofer-Lind safety index. 

The comparison of the C-L and RASCAL algorithms was made using the sample 
problem in Ref. 7. The problem is posed as the determination of the 
probability of failure of a cylindrical pressure vessel having an external 
torque. Failure occurs when the Von Mises stress exceeds the yield strength. 
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Figure 28 . Comparison Of QLM Model And Simulation Studies 


Thus, the failure function is defined as: 


G(u) = R - 


( 300* P 2 + 1.92*T 2 ) 0 - 5 


( 7 ) 


where 


u 

P 

T 

R 


Vector of random variables 

Pressure 
External torque 
Yield strength 


Each of the random variables is defined by the parameters in Table III. 
The vessel failure occurs when the value of G(u) is less than zero. 

The results of the RASCAL calculations are given in Table IV, and compared 
to the results presented in Ref. 7. As this Table clearly shows, the 
RASCAL method provides the same level of accuracy for the fa! lure 
probability calculation as does the C-L algorithm. In fact, the method is 
relatively insensitive to the RASCAL parametric values. There are several 
other important features of the calculations to be pointed out that 
indicate that the use of the RASCAL method is more appropriate for the low 
probability event calculations than is the C-L algorithm. The first of 
these is shown in Figure 29, where a plot of several of the RASCAL 
calculations is made versus the C-L calculation. The portion of the CDF 
shown in this figure shows the RASCAL calculations as single points for 
each individual run. Thus, for the RASCAL 250-40 where 250 is the number 
of intervals, and 250 times 40 (10,000) is the number of samples used, run 
each point denoted by an X represents the results of a single run of the 
RASCAL code. The straight line representing the C-L calculation is an 
interpolation (linear) between individual runs of the C-L program. This 
is necessary because the C-L algorithm provides only point estimates for 
the failure probabilities; it does not provide the entire CDF range, as is 
done in the RASCAL algorithm. To obtain the CDF of the failure 
probability it was necessary to run the C-L program 33 times. Of course 
this can be easily automated, however, it would still require a user to 
specify the levels at which the probability calculation be performed. 
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TABLE 9 

RANDOM VARIABLE DESCRIPTIONS 


VARIABLE 

DISTRIBUTION TYPE 

MEAN 

STANDARD 

DEVIATION 

R 

WEIBULL 

48.0 

3.0 

P 

LOGNORMAL 

0.9874* 

0.16* 

T 

EXTREME VALUE I 

20.0 

2.0 

‘Median 

and coefficient of variation 




TABLE 10 

COMPARISON OF FAILURE PROBABILITY CALCULATIONS 


METHOD 

PROBABILITY OF FAILURF 

MONTE CARLO 

1.600 x 10- 3 

CHEN-LIND 

1.820 x 10- 3 

RASCAL 10-50 

1.945 x 10- 3 

RASCAL 25-40 

1.819 x 10- 3 

RASCAL 20-400 

1.823 x 10- 3 


81 



At the end of the run, it .ay be discovered that a significant portion 
the probability curve was not covered, i.e. a portion of the COF 
missing. Again, this will lead to rerunning the program. 


Of 
i s 


In contrast the RASCAL method automatically takes care of covering the 
widest possible range of the CDF in the available computational time. or 
the RASCAL calculation in which the input variables were divided^into twen y 
discrete intervals, estimates of the CDF from the 10 leve of 
probability up to the 99.9999th percentile value is covered automatically by 
this algorithm. This range is covered only if the variables are independent. 


This study indicates 

probabilistic technique 


that the RASCAL algorithm is the most effective 

to use for generic space propulsion applications. 
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RASCAL Prediction: Failure Probability 


Exceedance Of Yield Strength 



Yield Strength — Von Mises Stress 

+ 10-100 o 40-25 X 250-40 Chin-Lind 


Figure 29. Rascal Versus Chen-Lind Failure Probability Predictions 




5.11 Comparison of ANLOAD Predictions With Expert Opinion. 


In the calculation of the engine loads during steady state operation, or 

when the loads are slowly varying, the probabilistic model contained in the 
ANLOAD computer program uses a set of engine influence coefficients that 
defines nominal operating conditions and effects of perturbations about the 
nominal point. The perturbations have been developed as a set of 33 
independent variables that are used to account for engine-to-engine and 
test-to-test variations on the SSME engine. 

As previously discussed, these variations were developed from consultations 
with the individual experts on specific hardware and covers geometric, 
performance, etc., conditions that can affect the engine operation. The 
purpose of developing these variations is to predict operating ranges for 
engine performance parameters such as pressures, temperatures, and flows. 
These random variations are added to the predicted performance effects of 

direct independent variation allowed by engine contract specifications to 
determine parameter minimum and maximum expected values. The direct 
independent variations include: propellant inlet temperatures and pressures, 
line resistance changes due to gimbaling, and tank repressurization flow 
settings. These maxima and minima define the operational limits used for 
engine component design and are used in developing the SSME engine balance. 

These types of engine variability estimates are the type of information 

required for generic load definition in the CLS code. The variability 
estimates combined with the engine test results furnish an ideal set of 
information for verification of this portion of the ANLOAD code. Only 23 of 
the 41 independent variables are used for the influence coefficient load 
calculation in ANLOAD. This limitation is not strict since a new set of 
influence coefficients, if supplied, could include up to 42 independent 
variables. This is only the method which is currently used for the ANLOAD 
program. Thus, the SSME balance calculate variations should be somewhat 
larger than those calculated with ANLOAD. Also, the engine balance 
variations are based on maximum operation limits, whereas the ANLOAD 

variability is based on actual engine test variability. This implies that 
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the engine balance calculations are more conservative than the ANLOAD 
calculations since they account for the design operating limits, while the 

ANLOAD results are only taking into account the variability actually seen 

during tests or standard operations. 

The first step in making the probabilistic predictions was to enter the 

expert opinion predictions of the variability in the independent parameters 
into the data base. These are reproduced in Table V. The coefficient of 
variation reported in Table V is the percentage value of the standard 

deviation of the nominal value at a specified power level (104%), assuming 
parameter in the engine, and should provide a reasonable test of the model. 

The variation listed in Table 11 is in addition to the variability induced 
by the random nature of the processes, e.g. engine to engine variation. 
Therefore, when the probabilistic function is constructed, it will have this 
variability as an additional parameter in the functional relationship. With 
the assumption that all of the independent parameters which control the 


85 



TABLE 11 

COEFFICIENT OF VARIATION FROM EXPERT OPINION 
FOR SSME INDEPENDENT VARIABLES 


Variable 


Coefficient 
of Variation 


Commanded Mixture Ratio 

HPFTP Turbine Efficiency Multiplier 

HPFTP Turbine Flow Multiplier 

HPOTP Turbine Efficiency Multiplier 

HPOTP Turbine Flow Multiplier 

T/C Characteristic Velocity Multiplier 

FPB Fuel Injector Resistance 

OPB Fuel Injector Resistance 

Oxidizer Pressurant Flow Rate 

LPFTP Inlet Orifice Resistance 

LPFTP Turbine Nozzle Area 

Fuel Pressurant Flow Rate 

LPOTP Pump Cavitation Correction 

HPOTP Pump Cavitation Correction 

LPFTP Pump Cavitation Correction 

HPFTP Pump Cavitation Correction 


0.5% 

1 . 0 % 
1 . 0 % 

1 . 0 % 

1 . 0 % 
0.125% 
1 . 0 % 

1 . 0 % 

1 . 505% 
1 . 0 % 

1 . 0 % 
0.65% 
0 . 0 % 
0 . 0 % 
0 . 0 % 
0 . 0 % 


86 


dependent variables are normally distributed, we can write, using the algebra 
of normal distributions and the influence function relationships: 


a m = E d m, i p 

(i-1) 


(8) 

s m “ E ( a m • 

b k . C0V k ) 2 

(9) 

“ E c k,m,i 

p(i-D 

(10) 


where d^. is the nominal dependent variable influence coefficients, 

c k,m, i is Uie 1n1 luer,ce C0<M i ic lent relating the independent variable k to 
the dependent variable m, C0V fc is the coefficient of variation of the 
independent variable k, a^ is the mean of the dependent variable m, and 
s m 1S the standard deviation of the dependent variable m. These equations 
were derived directly from the influence function equations and have been 
presented previously in the context of measurement error. Equations 8-10 
can be used to predict exactly the variability of any normally distributed 
set of independent parameters in the influence functions. 

The probabilistic model was run using these inputs to predict the 
variability in the HPFTP turbine speed. Table 12 gives the results of these 
calculations and Figures 30 and 31 present the results graphically. The 
first column gives the results when the mixture ratio was held constant at 
its nominal value. The second column of Table 12 gives the results when all 
of the variables in Table V were allowed to vary according to their 
respective coefficients of variation. 

These calculations were performed to show the validity of equations 8-10. 
If equation 9 is used to predict the variance when the coefficient of 
variation for the mixture ratio is set equal to zero then one finds that the 
standard deviation is equal to 45.6, which compares exceptionally well with 
the model calculation of 47.0 from Table 12. The plots in Figures 30 and 31 
show this difference, and they point out the dominance of the mixture ratio 
in determining the variability in the turbine speed. 
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originally projected. As expected the engine balance variation in e 
column of Table 12 do bound the test data results. 
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TABLE 12 

VARIABILITY IN THE HPFTP TURBINE SPEED 


PARAMETER 

MIXTURE RATIO 
HELD CONSTANT 

ALL VARIABLES 
RANDOM 

TEST DATA 

Mean value 

35500.0 

35502.5 

35425.3 

Computed standard 

47.0 

147.1 

341 .5 

Devi ation 




2- sigma 

0.26% 

0.83% 

1.93% 

Percentage 




Expert opinion 

1.10% 

1.10% 

1.10% 
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Figure 30. HPFTP Turbine Speed: Mixture Held Constant 
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Figure 31. HPFTP Turbine Speed: All Variables Random 



6.0 EXAMPLE TURBINE BLADE ANALYSIS 


6 . 1 Introduction And Definition?; 

To examine the possibilities in the application of the probabilistic model 
discussed in the previous section, four sample problems have been 
constructed. The first of these examines a change in the mission profile 
for the SSME HPFTP turbine torque analysis, from a steady state level of 
104% power to 109% power. The second example considers the changes in the 
SSME HPFTP turbine torque when the inlet turbine temperature is increased 
10% during the 109% power level operation. The third case is performed for 
tne SSML HPOiP turDine torque at 109% power. While it is true that this 
analysis also can be done based on previous data analysis, it will be 
performed using the probabilistic code, AN LOAD. In this way, the validity 
of such an approach can be examined. Finally, a prediction for the turbine 

torque for the turbine fuel pump in the J2 engine operating at 109% power 
will be examined. 

6 . 2 SSME HPFTP Tu rb i n e Torq ue At 109% Pnwer 

Table 13 gives the input variable definitions for this analysis. The 
procedure described in Appendix A was used to set-up the input to the 
probabilistic program for this analysis. The RASCAL analysis was used to 
calculate the HPFTP turbine torque with 1000 simulation points and each 
input discretized into 40 intervals. The results of the analysis Were a 
uniform distribution for the torque with a mean value of 10,824 ft- 1 b and 

a standard deviation of 64 ft-lb f . This indicates that the coefficient of 
variation is approximately 0.6%. 

6 • 3 SSMEjjP£T^TjJxbln, e ,.. Torque : 109% & 10% Increa se in Inl et Temperature 

For this study, the same inputs as given in Table 13 were used — except the 
mean value of the fuel inlet temperature was increased by 10%. It was 
assumed further that the coefficient of variation remained the same, so that 
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TABLE 13 

STANDARD INPUTS FOR PROBABILISTIC CALCULATIONS 


Vari abl e 

Commanded Mixture Ratio 
Fuel Inlet Total Pressure 
Oxidizer Inlet Total Pressure 
Fuel Inlet Temperature (R) 
Oxidizer Inlet Temperature (R) 
HPFTP Turbine Efficiency Multi 
HPFTP Turbine Flow Multiplier 
HPOTP Turbine Efficiency Multi 
HPOTP Turbine Flow Multiplier 
T/C Charac Velocity Multiplier 
Main Fuel Valve Resistance 
Main Oxidizer Valve Resistance 
Oxidizer Pressurant Flowrate 
FPB Fuel Injector Resistance 
OPB Fuel Injector Resistance 
LPFTP Inlet Orifice Resistance 
LPFTP Turbine Nozzle Area 
Fuel Pressurant Flowrate 


lyfie 

Mean 1 

Std Dev 2 

Uniform 

5.97443 

6.05108 

EV-I 

25.2313 

.173689 

Normal 

64.3341 

21 .0374 

Lognormal 

3.61308 

0.0162595 

Lognormal 

5.10174 

7 . 1 9274E-03 

Normal 

1.003 

.02018 

Normal 

1.0125 

.02025 

Normal 

1 .0152 

.020304 

Normal 

.9741 

.019482 

Normal 

1 .004 

.00251 

Normal 

.0138 

.0017526 

Normal 

.0107 

.0013589 

Normal 

0.0557 

0.0017267 

Normal 

.155 

.0031 

Normal 

.685 

.0137 

Normal 

.716 

.01432 

Normal 

.95 

.019 

Normal 

.032897 

.000428 


1 Thi s is the lower bound 
extreme value Type I. 

2 This is the upper bound 
extreme value Type I. 


for the uniform, and the 
for the uniform, and the 


shift parameter for the 
scale parameter for the 
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the standard deviation was calculated based on the coefficient of variation 
for the fuel turbine inlet temperature given in Table 13 times the new value 
of the mean. The results of this calculation indicate that the distribution 
is still uniform, and the standard deviation is still equal to 64 ft-lb 
but the mean value is now 10,937 ft-lb f . Thus, the assumption that the 
coefficient of variation for the fuel turbine torque is independent of power 
level, and therefore, is independent of time, is not valid. 


The reason for this can be seen easily by examining the model for the 
calculation of the torque. In the table of influence coefficients, the 
value of the fuel inlet temperature is not changed by simply changing the 
power level. Thus, while the change in the mean value will cause a change 
in the pieuictfeu torque value, it does not affect the variance. This is 
obvious from examining the equations for the influence function 
calculations. In equation (2b) for the fuel turbine inlet temperature, the 

c j,i ,k s are equal t0 2ero for k equal to 2, 3, and 4. Therefore, 
the variance contribution from changes in the fuel inlet temperature are 
equal to zero. 


6 . 4 SSME HPOTP_Torq ue Prediction From Scaling 

The derivation of a scaling parameter is based on the analysis just 
performed. Given the calculated value of the HPFTP torque, its horsepower, 
and its speed, the following scale parameter is used: 

T = 5300*HP*P L /RPM 

This gives a mean prediction of 5110 ft-lb f . Assuming a similar 
coefficient of variation gives a standard deviation of 30.7. 

Using the QLM to calculate the mean and standard deviation gives values of 
5083 and 35.0 ft-lb f , respectively. Obviously, these are relatively 
accurate. Additional studies must be performed to quantify this accuracy. 
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6.5 . 1 ? Fuel Turhine Tor q np Prediction From Scalin g 


Similar to the last case, the predicted value of the torque oh the 32 engine 
wild give a meah value of 1800 ft-lb f . If a value for the coef ticient : 
variation is used from the SSME HPFTP analysis, the predicted standa 
deviation is 11 ft-lb f . However, since the 03 is a single sta^e engine 
the variability from the inlet temperature should be reduce^ Si • 
contribution from the SSME HPFTP analysis is approximately 33% of the 
variance, it is predicted that the standard deviation should be e we 


and 11 ft-lbf. 
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7.0 LDEXPT : THE LOAD EXPERT SYSTEM FOR CLS 


7.1 Goal and Status 

The goal of the composite load spectra project is to provide a 
knowledge-based tool to generate and analyze composite loads of a rocket 
engine design and to supply them in a form that a probabilistic finite 
element computer program can use. 

This is being accomplished by developing probability models to simulate 
engine performance and other loads to collect the expertise built up over 
the years in order to help design a new and improved rocket engine. This 
computer program will provide a powerful probabilistic and statistical tool 
to guide users to obtain probabilistic information on rocket engine 
component loadings and provide expertise in analyzing engine loadings 
probabi 1 i sti cal ly . 

A knowledge-based system has the facility of building up a large domain 
knowledge base and maintaining a large amount of data. It has the 
capability to perform logical deduction and inferences and thus it can help 
users to make decisions and to solve problems. These characteristics allow 
one to build an expert system to simulate and perform the process of problem 
solving by an expert in a particular problem domain. 

This project requires a knowledge-based system that has a built in powerful 
probabilistic modeling and statistics tool box and a large database of 
rocket engine knowledge. This knowledge-based system will help the engineer 
to master probabilistic modeling technology and provide probabilistic 
information for structural analyses. The functions of this knowledge-based 
system are to manage the database, provide expert knowledge in generate 
probability loadings for rocket engine. In addition to being able to 
utilize the vast amount of existing FORTRAN probability and statistics 
tools, the code is required to have a FORTRAN based system. There is no 
existing knowledge system development tool that can satisfy all the needs of 
this project. Therefore, it was decided early that a FORTRAN based 
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non-proprietary knowledge system will be built to suit the needs of this 
project. 

A simple philosophy discovered by pioneering workers is that the power of a 
knowledge base system is in its capability to have a vast amount of domain 
knowledge and not necessary to have a complex inferencing engine. Following 
this philosophy, the load expert system LDEXPT was built with a simple 
inference system. This is the expert system driver controlling the rule 
processing and the user query interface. The expert system needs to know 
how to perform probabilistic modeling and statistical analysis. Therefore, 
a powerful probabilistic modeling and statistics tool box (consisting of 
FORTRAN routines) was built and is continuing to be developed. To make 
knowledge representation more efficient for the load expert system, a 
database system was implemented. This database system facilitates the 
communication between the expert system and the knowledge base, helps 
maintain data integrity and avoid data redundancy. The load expert system 
LDEXPT version 2.0 has all three elements in place. Its knowledge base has 
load information for SSME type engines, knowledge about the influence 
coefficient method for engine performance analysis and the turbine blade 
load information and scaling model calculation. 

The load expert system is a rule-based expert system. The inferences are 
carried out with the rules. In the load expert system the rules are 
modularized. Each module was designed to solve a particular problem or to 
perform a task. The load expert system LDEXPT version 2.0 has rule modules 
to calculate turbine blade loads using the scaling model and to generate 
engine dependent loads (e.g. HPFPT discharge pressure) using the influence 
coefficient method. The rules designed so far are mostly related to process 
control and information retrieval. In the next development, rules to 
generate probability models for a complicated composite load spectra will be 
designed which will require more intelligence. 

The probabilistic modeling and statistics tool box now has a stand alone 
load spectrum generator for steady and quasi-steady state loads. It has a 
statistics data analysis package which can select a best fit distribution 
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for a random variable and evaluate its distribution parameters. It also has 
a simple plotting routine to plot the duty-cycle-data profiles and a random 
walk, plotting routine to simulate a stochastic process. 

The knowledge base now has knowledge of the turbine blade loads for 
generating steady state and quasi-steady state load spectra. Additional 
load data on pressures and temperatures are ready for adding to it. The 
transient loads, pops and chugs and vibration loads are being developed and 
will be implemented as soon as the model development is complete. Knowledge 
on the transfer duct has been collected and rules for transfer duct load 
calculations can now be developed. 

The basic expert system components of the load expert system LDEXPT are all 
in place: the expert system driver, the database system, the FORTRAN data 
management system and the basic probabilistic modeling and statistics tool 
box; that is, the main tasks of system development phase are complete. The 
next main task is the implementation of additional applications of the 
expert system to the composite load spectra project. 

During the past year, the load expert system LDEXPT version 1.0 was 
implemented on IBM/PC and the NASA Lewis Research Center's mainframe. The 
IBM/PC version was implemented with Microsoft FORTRAN. The NASA/LeRC 
version was in IBM VS-FORTRAN. The IBM version has no database system and 
no plotting. When the LDEXPT version 2.0 is checked out on the Rocketdyne's 
Perkin-Elmer computer, it will then be implemented onto the NASA/LeRC 
system. This will have the database system and the plotting utilities. In 
addition, it will have additional rule modules for the turbine blade loads 
and the transfer duct loads. 
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7 . 2 IDEXPT. The Load Expert System 


A rule-based expert system casts knowledge into rules. It uses modus ponen 
and universal specialization to carry out the inference process. Rules are 
in the form of IF... THEN... which are called production. Thus, a rule-based 
system is also called a production system. In a pure production system, 
rules and control system (for searching) are in separate modules: the 

knowledge base and the inference engine. The knowledge base consists of the 
problem-solving knowledge, the process-control knowledge and the database 
knowledge, all in a rule-form. The function of the inference engine is to 
perform a chosen (built-in) searching algorithm on the knowledge (rule) base 
in order to reach a solution to one's problem. There are two types of 
searching strategies that are widely used in the rule base systems. They 
are the forward chaining and the backward chaining strategies. The forward 
chaining strategy starts from an initial problem state, searches forward 
until it reaches a goal state. The backward chaining strategy starts from a 
goal state, searches backward until it reaches the initial state. Searching 
forward means that rules are searched until the condition part of a rule 
(LHS, left hand side) matches with the present state, and then the 
conclusion part (RHS, right hand side) of the rule is used to move to a new 
state. Searching backward is just the opposite. Rules are searched until a 
rule is found such that the conclusion part of the rule matches the goal 
state. The condition part of the rule then becomes the new goal state 
(subgoal). This process is repeated until the initial state was reached. 

The pure production system is a very powerful tool. It is most suitable to 
the classification problems such as diagnostics. However, when problem 
becomes complex and deals with multiple data types at the same time, it is 
difficult to maintain an uniform knowledge base and the generic nature of 
the inference engine slows down the inference process. When frames were 
introduced, the generic inferencing scheme becomes very inefficient, if not 
impossible . 

The composite load spectra evaluation and generation is a planning and 
prediction problem. It has a large data base with multiple data structures 
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but is expected to have small number of rules for problem-solving knowledge 
and process-control knowledge. It needs to carry out numerous computations 
and it Is best to branch out of the expert system to do the analysis. It is 
obvious that a pure production system mechanism will not be able to satisfy 
the needs of this program. At the time, an expert system development tool, 
EXTRAN, was available to us for evaluation. It was a rule base system using 
a decision tree inference scheme. Its rules were built into a decision tree 
hierarchy. Searching through the tree was carried out with user supplied 
information or selection. The tool was in FORTRAN, which had the advantage 
of a convenient integration with vast amount of the probabilistic modeling 
and statistics tools available for engineering. The loss of flexibility of 
a pure production system was not serious because the rule base for the 
problem-solving knowledge and the process-control knowledge of this program 
was expected to be not too large and rules would not be modified frequently 
once they were established. So the load expert system inference scheme was 
conceived and it was modeled after EXTRAN. There is an added benefit to 
have the load expert system compatible with EXTRAN. EXTRAN is used at 
Rocketdyne to develop expert systems for engine performance analysis and 
ig frequency data analysis. These different programs could benefit each 
other by using similar expert system development tools. 

The load expert system for the composite load spectra project, named LOEXPT, 
has the rules built in a decision tree routine or several routines if the 
rules can be decoupled. In this way, the rules are modularized. The load 
expert system communicates with users via a problem text files where 
questions for query are stored. The system records the process so that it 
can show them to the user when a "HOW" is asked. User can also ask “WHY" a 

question is prompted and the load expert system will reveal the logics 
behind the question. 


The load expert system LOEXPT version 1.0 was built with the scheme 
escribed above. A very simple consultation system was built and 
demonstrated the feasibility of the design. However, the system needed a 
ot of help from users in the following ways. Many redundant queries had to 
be made to obtain information which could be identified with one or two 
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Figure 32. LDEXPT: Load Expert System 
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7.3 Implementation of the Load Da t abase System and Interface 

The engineering data in LDEXPT version 1.0 was organized in a conventional 
data file format. This format was found to cause data redundancy problems, 
which the expert system could not handle. For example, if a test-data ID 
number was identified, the engine type and mission history profile type were 
determined (known to our engineers). However, the expert system would 
require additional specific rules to identify these relations. This 
redundancy problem not only complicated the system but also resulted in an 
exponential growth in the number of rules required for the expert system. 
This problem could be resolved if the engine type and the mission type were 
built into a property list of the test-data ID as could be done in LISP or 
if the three attributes were built into a database, e.g. a relational 
database. Using test-data ID as a key, once it was identified, the engine 
type and the rest of the properties were determined. The expert system did 
not need to query further to acquire other information. 

Organizing data into a database model has many advantages. The obvious ones 
are avoidance of data redundancy and inconsistency, ease of enforcing data 
integrity and ease of data maintenance. A database model of engine data 
also has the side benefit of a well organized data base and an easy to 
understand retrieval system. Moreover, the most important advantage of 
building a database is that it facilitates the communication between the 
expert system and the engineering data base, which in turn speeds up the 
knowledge acquisition process for the load expert system. 

A Relational Database Model . A relational database model is like a table of 
key variables and attributes. The keys are used to identify a record (row) 
in the database model (table). The values of the keys are unique to each 
record. In database terminology, the model is in normal form. For the load 
database, for example, a- database table LOAD could be built which has the 
following fields (columns): 
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LOAD-ID, LOAD-name, Mean, COV, Di st-type, NE-coeffs 


where LOAD-ID is the key, 

LOAD-ID is the load ID number 

LOAD-name is the load variable name, e.g. mixture ratio 
Mean is the default mean value of the load at 100% power 
COV is the coefficient of variation of the load 
Dist-type is the default distribution type of the load 
NE-coeffs is the inference coefficients for calculating 
the nominal engine mean value of the load at the desired 
power level 

With this database table, once the load ID is identified, the mean of the 
load, the COV and the rest of information can be easily accessed by the 
expert system. 

A second example is the Duty-Cycle-Data table: 

Test/FI ight-ID, Load-ID, Engine-type, Mission-type, Duty-Cycle-Data 

where Test/FI ight-ID and Load-ID are the keys 

Test/FI i ght-ID is the test data ID or flight data ID 
Load-ID is the independent load ID 
Engine-type is the engine type, e.g. SSME 
Mission-type is the mission history profile type, e.g. 
acceptance test 

Duty-Cycle-Data is the group-name of the duty cycle data 
stored in the data file 

There are a total of five databases built for calculating the turbine blade 
component load spectra: LDIP, LDEP, INFC, LTBC and DFAT. They are listed 

in tables X.l to X.5 of Appendix B (for INFC, only samples are listed 
because of its large size). The database description of the five tables are 
also listed in Appendix B. 
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The three databases LIDP, LDEP and LTBC belong to the same class of the 
object LOAD, which possesses the following attributes: load-ID, load name, 
mean, coefficient of variation (COV), P3 (rare event probability limit), 

distribution-type, and a set of coefficients for its nominal engine 

configuration. In the load databases, the mean values for the loads are 

values for the nominal engine configuration. The COV's and their 

distribution type for mixture ratio, fuel inlet pressure and temperature, 

and LOX inlet pressure and temperature are values based on engine data. The 

COV's for the rest of the independent loads and all of the dependent loads 
are values based on expert opinion and the SSME engine balance model. The 
distribution type for all other loads was assumed to be normal. The nominal 
engine coefficient sets were obtained from the influence coefficient file 
"INFLUENCE.DAT". These entries, that is the mean, the COV etc., for the 
turbine blade component load are not available at this time. When default 
values are available, they will be stored into the LTBC database. 

The INFC database has sixteen (16) tables. Each table includes information 
for four dependent loads. For example, the first table has influence 
coefficient set and gains for dependent load 1 to 4 , the second table is for 
dependent loads 5 to 8 and so on. an independent load. The attributes for 

INFC are the dependent load-ID, the independent load-ID, the influence 

coefficient set (4 numbers) and the gain set (4 values). The influence 
coefficient sets were obtained from the influence coefficient file. The 
gain set includes GAIN65 (gain for 65% power level), GAIN90 (gain for 90% 
power level), GAIN100 (gain for 100% power level) and GAIN104 (gain for 104% 
power level). These gains were calculated based on the assumption of normal 
distribution for the independent load and using the COV values from the 
independent load database LIDP. The idea of including a gain set covering 
the operating range of power level is significant. By examining the values 

of the gain set, one could easily spot which dependent load gain was 

significantly contributed to by the independent load. This kind of expert 
knowledge can now be built into the load expert system thanks to the 
database system. Here we have learned an important lesson that is: the 

knowledge representation of the domain knowledge is very important to the 

success of this project. 
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Database Design . A relational database model is being built for LDEXPT. An 
indexed sequential access method (ISAM) algorithm is employed for retrieval 
of database records. A key file is constructed for each database table. The 
keys are sorted in certain order. The records are then retrieved through 
the index stored in the key file uniquely identified by the values of the 
keys. There are physically two files for each table, a data file contains 
all records of data and a key file contains values of all the keys. In this 
model, no secondary keys are allowed. The different key variables in the key 
file are variables of a primary split key, all values of the split key 
variables must be identified uniquely for a record retrieval. 

Functions identified fo r the database system are CRFATF a table, UPDATE a 

table, DELETE a record, SELECT a record, BUILD a key file and SAVE a table. 
The detail function of each procedure is presented below. 


CREATE a table: 

Set up data dictionary: record description: # of field, # of keys; 

field description: field-name, data-type; 

Enter data records; 

BUILD a key file; 

Save the database. 


READ data dictionary: 

Open a database file and locate the desired database table; 

Read data dictionary (i.e. field-names or key-names of the table); 
Move database table data on-line. 


UPDATE a table: 

Enter a new record; 

Sort the new record key(s) into the key file; 
Save the database (optional). 
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DELETE a record: 

Delete the record key value(s) from the key files; 
Mark the delete data record in the table. 


SELECT record(s) : 

Retrieve data record index/indices from the key files; 
Di splay the record(s) . 


BUILD a key file: 

Create a key file; 

Sort the key value(s) . 


SAVE a database table: 
Select SAVE option; 

Save the database table. 


A very efficient algorithm for building an ISAM key file is B-tree. 
However, because there is no pointer data type and no recursion in FORTRAN, 
it is difficult to build a B-tree, in fact any tree. For now, a sorted 
array will be used as our key file format. 

The database functions being built for LDEXPT are similar to those 
identified above. The available database commands and the database routine 
descriptions are presented in Appendix B. 


Database Limitation . For the moment, the database model is limited to 15 
fields, 10 (split) keys and 100 records. Each field of character type is 
limited to 8 characters long. 
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Keys are only of character or integer type. These limits were chosen based 
on the needs of LDEXPT and they can be expanded easily if they do not become 
too large. For building large size database, this model has to be modified 
and the limitation of using FORTRAN language to build a database will 
severely hinder the effort. Other languages such as C would be more suitable 
for that purpose. 

The main purpose of this database system is to enable one to write rules to 
retrieve knowledge for the expert system. Hopefully this knowledge can be 
modeled with many small relational database tables. Normalized relational 
databases tend to be comprised of small tables. Knowledge in this form will 

be easier to be understood by engineers who have to sort out large volume of 

data. 

Interface . Communication between the expert system driver and the database 
system is achieved by putting an expert system option into the database 
routines. The interface allows the expert system driver to query only the 
key attributes and to retrieve data items from a database. Two interface 
routines GRSPRC and GRSPMN were written. GRSPRC grasps (retrieves) a 
database record and stores it in the array ITEMS. GRSPMN grasps an item in 
a database record. These two generic data retrieval routines send the 
desired data items to the rule base system for processing. 

The implementation of the database system is a significant development for 
the load expert system LDEXPT. With appropriate representation, the text 
book knowledge and expert's knowledge can be incorporated in the knowledge 
base to make LDEXPT into an intelligent load expert system. 
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7.4 LDEXPT's Rules and Implementation 


The domain knowledge for the composite load spectra (CLS) project consists 
of two main areas: the probabilistic modeling method and the rocket engine 

structural load information and calculation. The synergism of the two 
domains have to be brought about to produce the domain knowledge for CLS. 

Knowledge acquisition is the key for building a successful expert system. 
The rocket engine domain knowledge covers a broad range of information: 
rocket engine component geometric information and operating condition, 
rocket engine measurement such as engine performance data and power spectral 
distribution, rocket engine structural load models etc. There is a rich 
pool of information built up from the last several decades. Some are in 
notebook form, others are in textbooks, and many of them are measured data 
stored in data files, in LOTUS files or simulated in models such as the 
engine balance model and the influence coefficient model. There is also a 
vast amount of knowledge built up over the years in the rocket engine 
specialists minds. Many of these expertises are not documented anywhere. 
This knowledge is being derived from specialists at Rocketdyne. We are 
consulting with specialists who work on the on-going SSME data collection 
and evaluation tasks to supply data information relevant to this project. 
We find out how experts use models to simulate engine performance. All 
these have to be built into an uniform framework so that the expert system 
can utilize them effectively. The framework for storing knowledge is 
another key for building a successful expert system. 

The probabilistic modeling is the other important knowledge domain for the 
CLS project. It includes modeling random variables, simulating stochastic 
processes, data statistical evaluation and many more. Battel 1 e ' s expertise 
in this area is being used to help build models to simulate composite load 
spectra. Most of these 'knowledges are in algorithmic forms and could be 
built into computational procedures. The remaining problem is of course 
that how does the expert system communicate and utilize these procedures. 
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All the knowledges relevant to the load expert system come in three basic 
types: text information, engineering data and modeling algorithm. Rules 

are designed to represent the textual information needed for load 
generations. Rules also are written to control the computation process and 
to retrieve and manage the requested engineering data. For the load expert 
system LDEXPT, rules are separated into rule modules to perform specific 
tasks. Interactions between different rule modules are controlled by other 
modules which employ a simple working memory model to communicate between 
them. The model will be improved throughout the project. 

The Working Memory Model . The working memory model was designed for passing 
information (short-term memory) between different rule modules. To keep the 
model simple, the information saved was limited to that needed to pass from 
one module to another module but not between multiple rule modules. The 

working memory consists of a "stack" and a memory array. The "stack" is 
used for storing database indices for record retrieval and the memory array 
is used for storing information (e.g. subgoals, facts). 

The advantage of implementing a working memory model is that many inference 
processes can proceed without user intervention. For example, suppose a 
turbine blade HPFTP centrifugal load spectrum calculation was requested. 
The load expert system would first find out what dependent load(s) was 
needed and the associated scaling model coefficient from the turbine blade 
load database. In this case, it was the HPFTP turbine speed. This 

information was stored in the working memory and passed to the dependent 
load rule module. There the expert system retrieved the dependent load 

information and the associated nominal engine coefficient sets. Then, the 
expert system began to select the independent loads with the help of another 
rule module. The dependent load ID's were passed from the dependent load 

rule module. The expert system then selected a number (as requested by 

user) of independent loads which had the most contribution to the turbine 
speed load. Without the working memory, user has to supply this information 
between modules in order to complete the process. 
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mFXPT's Rules . The rule base for the load expert system LDEXPT was 
redesigned since the implementation of the database system. The rule base 
is composed of rule modules. Each rule module is an independent unit and 
has its own database table(s) associated with it. Using rule modules as 
building blocks allows incremental development of a rule base for the expert 

system. 

The following examples are two of the rule modules implemented. 

Rule module for load data base: 

If the lead ID is number N 

then its name is AAAA, 
its mean is M 

its standard deviation is SIGMA, 

its P3 (rare event probability) is ZERO or a fraction, 
its nominal engine coefficients are A1 , A2, A3, and A4. 

Rule module for Influence Coefficient Model: 

a) If the dependent load ID is M and 

the independent load ID is N 

then the influence coefficient parameter set is Cl, C2, C3 and C4. 

b) If the dependent load ID is N and 

the user requests that the expert system selects the M most 
influential independent loads on the dependent load and 
the selection is going to be based gain for power level X 
then the M independent loads are Ml, M2.... 

c) If the dependent load ID is N and 

the user requests a simple deterministic influence 
coefficient model calculation 

then the expert system will either request user to select the 

independent loads manually or the expert system will select them, 
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retrieve influence coefficient set and 

perform the deterministic influence coefficient model 

calculation. 

Nine rule modules were implemented. These include rule modules for the load 
database, the duty-cycle-base, the influence coefficient model, the quick 
look model and the turbine blade load scaling model. Their rules are listed 
in Appendix C. 

More rule modules will be designed and built. It can be seen that the load 
expert system now possesses good knowledge about rocket engine loads of 
which some data (e.g. most of the load's standard deviations) are based on 
the expert's estimates. The load expert system also possesses good 
knowledge on the influence coefficient model and the turbine blade component 
load scaling model and will include the information on the other three 
components as the system progresses. It even shows signs of intelligence in 
being able to select the most influential independent loads for a dependent 
load calculation. 

There is no new information supplied to the expert system to enable it to 
select independent loads for the users. The information is in the data file 
INFLUENC.DAT (the data file used to perform influence coefficient model 
calculations). However, if one looks at the data file, one sees lines and 
lines of numbers. It is very difficult to extract any information out of 
it. After processing it and putting it into a load database, the 

information becomes alive. Simply by sorting the gain values that were 
built into the database, the expert system is able to select the most 
influential independent loads for a dependent load calculation. The point 
is that an appropriate knowledge representation is very powerful. The 

database implementation to the load expert system proves to be a very 
significant contribution to this development. The database and the rule 

modules from our load expert system together look very much like the frames 
used in the frame-based knowledge system. However the object programming 

paradigms such as inheritance and class so prominent in a frame-based system 
are not seen in our system. 
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The rules designed so far are mostly related to process control and 

information retrieval. As we acquire more and more expertise on how an 
expert solves a composite load spectra problem, we will build rules to carry 
out the induction process of the experts and in so doing increase the 
intelligence of the load expert system. 

7.5 l DEXPT Future Development 

The new version of the load expert system LDEXPT (version 2.0) is not yet 
fully tested. Different consultation sessions will be run to verify the 
correctness of the rule modules implemented. In the coming months, 

transient model and nonstationary stochastic algorithm will be implemented. 
More rules for turbine blades, transfer ducts, lox posts and the HPOTDD will 

be added. The pops, chugs and vibration data will be transformed into 

databases and stored in the load knowledge base. Probabilistic models for 

generating pops, chugs and vibration loads will be provided in the load 

expert system. 

The basic elements of the load expert system are all in place. The main 
task now is knowledge engineering. This involves designing representation 
for the vast amount of engine data, implementing the process-control 
knowledge and learning the problem-solving knowledge from experts and 
translating it into rules. As mention earlier, the two major knowledge 
domains for the composite load spectra problem are the rocket engine 
analysis and design, and the probabilistic modeling of loads. The experts 
from Rocketdyne and Battelle are relied upon to acquire the knowledge. 
Engine data and load information are analyzed and cast into a convenient 
form such that the load expert system can utilize them to perform 
intelligent tasks. Rules for the problem-solving knowledge need to be 

identified and implemented. 

The way the load expert system was set up requires that a knowledge engineer 
who is familiar with the system maintains the expert system. New knowledge 
and learning are added to the system manually by the knowledge engineer. 
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The procedure is not difficult. It needs to put the new data into a 
database and to write new rule module (FORTRAN) routines, compile them and 
link, into the load expert system. However, the knowledge engineer has to 
check the consistency of the new database with the existing ones, avoid 
implementing redundant data. He should also beware of any conflicts between 
rule modules and try to resolve them. Only then a sound expert system is 
maintain as it grows larger and more intelligent. 
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APPENDIX A 


PROBABILISTIC MODEL DRIVEN CODE FOR STAND ALONE OPERATION 


There are two primary computer programs for executing the probabilistic load 
model: (1) a BASIC program (MENU) for assisting in generating input and 
displaying the results of the probabilistic load; and (2) a FORTRAN program 
(ANLOAD) for performing the actual probabilistic load calculations. The 
steps for executing these programs are discussed below 

A program, MENU, has been written in BASICA that performs two functions. 
First it is used to generate an input file for use with the ANLOAD program. 
This is a menu driven program that writes an output file to a floppy disk 
which is subsequently used as input to ANLOAD. The second function of MENU 
is to display the results of the ANLOAD calculations. 

To illustrate the use of these programs, the sample problem used is the one 
discussed in last year's annual report. In this example there are six 
stochastic loads which are all stationary. The object is to calculate the 
mass flow rate in the HPFTP and the HPFTP shaft speed. Each of these 
quantities are combinations of four other, individual loads. 

To begin the session the IBM is turned on and the following commands are 

entered. It is assumed throughout the discussion that at least two disk 
drives are avai lable. 

C> GRAPHICS 
C> BASICA NEWMENU 

The MENU program is now operating. To use the MENU there are two important 
keys - the function key FI and the cursor down arrow. Simply stated the FI 
key can be thought of as a negative response (reject the option) while the 
cursor down key is a positive response (accept the option). The first 
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message on the screen is a "billboard" message. A billboard message does 
not require a response but rather is simply providing information to e 
user. Striking the FI key will cause the next menu to printed. 

The first message informs the user that the program does perform two 
different functions; either input preparation or display the lo» d mo e 
calculational results. After striking the FI key the program asks that 
of these two functions be selected. Since we are interested in preparing 
ioad model, run the first option is selected by striking the cursor down 

arrow. 

The nevr message lists all of the possible engine parameters and individual 
loads currently programmed. To include one of these variables in e 
calculation, the cursor down key is used; to not include one the 1 functio 
key is used. When all variables have been selected the F3 function key is 

pressed . 

The next screen lists the selection and asks if the list is correct. If it 
is the cursor down key is used, if not the FI key will return the user to 

the previous menu. 

The next menu asks which of the three probabilistic models is to be used for 

, j 4 if h <<• desired to use medium or high 

the current calculation and, then, if it is aesireu 

accuracy. 

The next set of menus allows the user to define the probabilistic form for 
each of the input variables. By pressing the F3 function key first the user 
requests that default variables be used. Otherwise, the cursor down an FI 
function keys work as usual. The program then lists the input and allows 
the user to accept or reject it before proceeding to the next variable. 

At this point all of the necessary data for the Independent variables has 
been selected. The next menu lets the user select the dependent variab es 
!h-ch he wishes to include in the analysis. The current dependent variables 
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are listed and the cursor down and FI keys are used to select or reject the 
variables on the list. The F3 key lets the program know that the selection 
process is finished. After the F3 key is used, the program lists the 
dependent variables to be included in the analysis and allows corrections to 
be made if there is an error in the list. 

Next, for each dependent variable the list of independent variables included 
in this analysis is presented on the menu. The cursor down key is used to 
select the variables from this list which are variable for the specific 
dependent variable being shown. The FI key deletes that independent 
variable from the input to the current dependent variable, and that variable 
only. If the independent variable is not to be included for subsequent 
analyses it must be deleted later. 

The next menu asks if there are any rarely occurring loads which would be 
added as "spike" type loads. If there are, the frequency of such loads is 
requested as well as what its form is. 

Finally the program requests information on the mission phase history. For 
the initial phase the type (transient, quasi-steady, or steady), its start 
time and the power level at this time, and its end time and power level are 
requested. Subsequent requests will not request the start time since it is 
assumed to be a continuous process in time. 

After all of this information has been collected, it is stored in a file 
named ANLOAD.DAT. This file can then be input directly to the FORTRAN 
program for analysis. 

A sample input is shown in Table A-l . 
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Table A-l 


Sample Input To ANLOAD 


5 

1 

2 

4 

12 3 4 

C 

100 1 5 999999 

2 

50 50 

HPFTP Turbine Speed 
Commanded Mixture Ratio 

1 5.97443 6.05108 0 

Fuel Inlet Total Pressure 

2 28.5545 7.38417 0.001 

Oxidizer Inlet Total Pressure 

2 64.3341 21.0374 0.001 

Fuel Inlet Temperature (R) 

3 3.61308 .0162595 0.001 

Oxidizer Inlet Temperature (R) 
3 5.10174 7 . 1 9274E-03 0.001 


4 3 0 0 

1 0 .65 

2 39000 500 

2 2.5 .1 

2 .65 1.04 

3 1.04 1.04 

100 


0 .05 

.05 5 


5 20 
20 100 
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APPENDIX B 


THE LDEXPT LOAD DATABASE DESCRIPTION AND EXAMPLES 


LIDP : independent load database table group name 


FIELD NAME 

DESCRIPTION 


LIDP-ID (key) 

independent load ID 


LD-NAME 

load name 


MEAN 

nominal engine mean value of the 

load 


or any predetermined default mean 

val ue 

cov 

default coefficient of variation 


P3 

rare event probability limit 


DIST 

distribution type 


NE-COEF1 

nominal engine coefficient, Oth order 

N E-COE F2 

1st order 


N E-COE F3 

2nd order 


NE-C0EF4 

3rd order 


LDEP : dependent 

load database table group name 


FIELD NAME 

DESCRIPTION 


LDEP-ID (key) 

dependent load ID 


LD-NAME 

load name 


MEAN 

nominal engine mean value 


COV 

default coefficient of variation 


P3 

rare event probability limit 


DIST 

distribution type 


NE-COEF1 

nominal engine coefficient 


N E-COE F2 

same as above 


NE-COEF3 

same as above 


NE-COEF4 

same as above 
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INFC : influence coefficient database table group name 


FIELD NAME 

DESCRIPTION 

LDEP-ID (key) 

dependent load ID 

LIDP-ID (key) 

independent load ID 

INFL-C1 

Oth order coefficient of the influence 
coefficient set 

INFL-C2 

1st order coefficient 

INFL-C3 

2nd order coefficient 

INFL-C4 

3rd order coefficient 

GAIN55 

unit gain for power level at 65% 

GAIN90 

unit gain for power level at 90% 

GAIN 100 

unit gain for power level at 100% 

GAIN 104 

unit gain for power level at 104% 

LTBC : turbine blade 

component load database group name 

FIELD NAME 

DESCRIPTION 

TB-C-ID (key) 

turbine blade component ID 

TB-LD-ID (key) 

turbine blade component load ID 

TB-LD-NA 

turbine blade component load name 

MEAN 

default mean value 

COV 

default coefficient of variation 

P3 

rare event probability limit 

LD-TYPE 

T/B component load type, e.g. point, 
distributed or nodes 

LDEP1-ID 

dependent load for scaling model calculation 

LDEP2-ID 

dependent load for scaling model calculation 

SC-COEF 

scaling model coefficient 

set to zero if more than one coefficients 

required 

TBC-GRPN 

T/B scaling coefficient file group name 
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DFAT : SSME flight and test duty-cycle-data file group name 


FIELD NAME 

DESCRIPTION 

DCD-ID (key) 

SSME flight or test data ID 

LIDP-ID 

independent load ID 

set to zero for engine power duty-cycle-data 

ENGINE 

engine type 

MISSION 

mission history profile type 
e.g. flight or acceptant test etc. 

DCD-GRPN 

duty-cycle-data file group name 
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TABLE B . 1 


LI DP: INDEPENDENT LOAD DATABASE 
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TABLE B.2 


LDEP: DEPENDENT LOAD DATABASE 
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TABLE B.2 (CONTINUED) 
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TABLE B. 3 


INFC: INFLUENCE COEFFICIENTS AND 

GAINS DATABASE (SAMPLE, GROUP #1) 
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TABLE B. 3 (CONTINUED) 
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TABLE B . 3 ( CONTINUED) 
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T8-C-10 


: TABLE B.4 

LTBC: TURBINE BLADE COMPONENT LOAD DATABASE 


to* ID 

T8-LD-NA 

LD-TYPE 

IDEPI-ID 

IDEP2* ID 

SC- COE F TBC-C 

1 

HF-CFO 

POINT 

2 

0 

0.10000E+01SCF1 

2 

H F - HP • T 1 

POINT 

46 

0 

0.11623E-01SCF1 

3 

HF-KP-T2 

POINT 

46 

0 

0.19194E-01SCF1 

4 

KF-HP-A1 

POINT 

54 

62 

0.98384E-01SCF1 

5 

HF-HP-A2 

POINT 

54 

62 

0.78707E-01SCF1 

6 

HF-T-11 

OIST 

46 

0 

0.11623E-01SCF1 

7 

HF-T-T2 

D1ST 

46 

0 

0.19194E-01SCF1 

ft 

HF-T-A1 

DIST 

54 

62 

0.98384E-01SCF1 

9 

HF-T-A2 

D1ST 

54 

62 

0.78707E-01SCF1 

10 

HF-HW-T1 

DIST 

46 

0 

0.11623E-01SCF1 

11 

HF-MM-T2 

OIST 

46 

0 

0.19194E-01SCF1 

12 

HF-HN-A1 

OIST 

54 

62 

0.98384E-01SCF1 

13 

HF-HH-A2 

DIST 

54 

62 

0.78707E-0;sCFt 

14 

HF-H-T1 

OIST 

46 

0 

0.11623E-01SCF1 

15 

HF-H-T2 

DIST 

46 

0 

0.19194E-01SCF1 

16 

HF-H-A1 

OIST 

54 

62 

0.98384E-01SCF1 

17 

HF-H-A2 

OIST 

54 

62 

0.78707E-01SCF1 

18 

HF-T-DP 

NODES 

19 

46 

0.00000E+00SCF1 

19 

HF-HN-DP 

NOOES 

19 

46 

O.OOOOOE+OOSCF1 

20 

HF-H-DP 

NODES 

19 

46 

O.OOOOOE+OOSCF1 

21 

HF-DYN-P 

POINT 

2 

0 

O.OOOOOE+OOSCF1 
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TABLE B.5 


DFAT: FLIGHT AMD TEST 

DUTY-CYCLE-DATA DATABASE 


DCO-ID 
STS61-A 
STS61-A 
STS61-A 
STS61-A 
STS61-A 
STS61-A 
902-38; 
902-38; 
902-384 
902-384 
902-3a*i 
902-384 
902-387 
902-387 
902-387 
902-387 
902-387 
902-387 
750-262 
750-262 
Hit any key 


LIDP- ID ENGINE 

0 SSME 

1 SSHE 

2 SSHE 

3 SSHE 

4 SSHE 

3 SSHE 

0 SSHE 

1 SSHE 

2 SSHE 

3 SSHE 

4 SSHE 

5 SSHE 

0 SSHE 

1 SSHE 

2 SSHE 

3 SSHE 

4 SSHE 

5 SSHE 

0 SSHE 

1 SSHE 
to continue 


> 


MISSION 

DCD • GRPN 

FLIGHT 

PURI 

FLIGHT 

KXR1 

FLIGHT 

PI FI 

FLIGHT 

PI01 

FLIGHT 

T I FI 

FLIGHT 

TI01 

TEST 

PUR 3 

TEST 

MXR3 

TEST 

PIF3 

TEST 

PIQ3 

TEST 

TIfi 

TEST 

TJ03 

TEST 

PWR4 

TEST 

KXR4 

TEST 

PIF4 

TEST 

P104 

TEST 

TIF4 

TEST 

TICK 

TEST 

PWR5 

TEST 

MXR5 


DCO-ID 

750-262 

750-262 

750-262 

750-262 

901-495 

901-495 

901-495 

901-495 

901-495 

901-495 

901-491 

901-491 

901-491 

901-491 

901-491 

901-491 


HOP-ID 

2 

3 

4 

5 
0 
1 
2 

3 

4 

5 
0 
1 
2 

3 

4 

5 


ENGINE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 

SSHE 


MISSION 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 


OCD-GRPM 

PIF5 

PI05 

TIF5 

TI05 

PWR6 

MXR6 

PIF6 

PI06 

TIF6 

TI06 

PUR 7 

MXR7 

PIF7 

PI07 

TIF7 

TI07 


129 



APPENDIX C 


database commands and routines description 


Available Patabas_e_CoiMands: 


7DBCR : 
7DBBK : 
7DBDF : 
7DBSL : 
7DBDL : 
7DBUP : 
7DBRD : 
7DBSV : 
7DBLT 
7DBLK 
7DBCF 
7INLD 
7INFL 


Qfeats a database table 
Build a database key file 
Display field and key names 
Select database records 
Delete database records 
Update (add) database records 

Open a database file and read its data dictionary 
Save a updated database table 
List a complete database table 
: List a complete database key file 
: Create fields for a database table 

. i I — J 


Build a load ID and properties database 

r • _ * i. _ nai 


d.,41 a influence coefficients and 




DBMS Routines: 


DBMS 

DBCRTB 

DBBUKE 

DBDPkF 

nRRDDC 


Database System driver 

7DBCR command, create a database table 

7DBBK command, building a key file 

">DBDF command, display field and key names 

inRRn command, open a Database file and read a Database 


table dictionary 
DBSLRC : 7DBSL command, select, 
Database table 


retrieve and display records of a 
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DBUPTB : ?DBUP command, update (add) new records to a Database table 

DBDLRC : ?DBDL command, delete records from a Database table 

DBSVTB : ?DBSV command, save a Database table to a Database file 

DBLSTB : ?DBLT command, list data on a Database table 

DBLSKF : ?DBLK command, list key data of a Database table 

DBGEIN : get record indices for the requested records 

DBWRFD : display selected records on CRT 

DBWRRC : write a retrieved record to CRT 

PRPAGE : print a page of data on CRT 

DBRDKD : read field and key descriptions from terminal input 

DBRDDA : read Database table input from terminal 

DBRDFD : read a field data from terminal input 

SORKEY : set up a multiple sort procedure and call SHLS02 

SHLS02 : a shell sort routine for a two-column-array 

DBSWIT : switch (substitute) row(indexl) by row(index2) 
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APPENDIX D 


LDEXPT Rule Modules and Routine Descriptions 


Rule module for duty-cycle-data base: 

If the flight or test data ID is XXXXXXX and 
the independent load ID is N 
then the engine type for this data is YYYY, 

thp mission history profile type is ZZZZZ, 

the data is stored in the duty-cycle-data base file with 

group name AAAA. 

Rule module for the Quick Look Model: 

If the dependent load ID is N and 

the user requests a quick look model calculation 
then the expert system will either request user to select the 

independent loads manually or the expert system will select 

them, 

retrieve influence coefficient set and 
perform the quick look model calculation. 

Rule module for turbine blade load scaling model 

a) If the turbine blade component ID of interest is M and 
the turbine blade component load ID is N 
then the turbine blade component load name is AAAAAAA, 
its load type is.TYPE-X, 

the dependent loads needed for the scaling model are load ID1 

and load ID2 (if required), 

the scaling model coefficient is NNNN or 

the coefficients are stored in a duty-cycle-data base file 

with group name BBBB (if more than one coefficient is 

required) . 


132 



If a turbine blade load calculation is requested and 
the turbine blade component and load are M and N 
then the expert system will generate an input file for an ANLOAD 
calculation which includes the following information: 
the independent and dependent loads required and 
their relevant load parameters and 
the influence coefficient sets and 

the duty-cycle-data for dependent loads if necessary and 
computational parameters such as number of bins, time slices 
etc. by prompting the user. 

t *p fnrhi no M r o m p n ne n t snd lo^d m ?. n d N 

the user requests a simple turbine blade scaling model 
calculation 

then the expert system will retrieve the default dependent load 
information and scaling coeffi ci ent(s) and 
perform a turbine blade scaling model calculation. 



The rule module routines: 


RBIDPL : rule module for retrieving independent load 

information and selecting independent loads for users 
based on the gain database 

RBDEPL : retrieving dependent load information manually or by 
the expert system with the help of the simple working 
memory model 

RBTBCL : retrieving turbine blade component load information 
and scaling model information 

RBQLM : the quick look model, calculating dependent loads 
assuming all loads are normally distributed 

RBSICM : the deterministic influence coefficient model, 

calculating point values for dependent loads using 
influence coefficients 

RBDRIV : the new rule base driver routine 

RBICGN : retrieving the influence coefficient set and the gain 
database 

RBSSM : rule module for performing simple scaling model 
calculation using default dependent load values 

RBDCD : retrieving flight or test duty-cycle-data files 

RBTBIN : preparing an AN LOAD input file for a full blown AN LOAD 
calculation 
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